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ABSTRACT 

A study of outflow (OF) and BAL systems in Mrk231 (and in similar IR QSOs) is 
presented. This study is based mainly on ID and 2D spectroscopy (obtained at La 
Palma/WHT, HST, IUE, ESO/NTT, KPNO, APO and CASLEO observatories) plus 
HST images. For Mrk 231 we report evidence that the extreme nuclear OF process has 
at least 3 main components on different scales, which are probably associated with: 
(i) the radio jet, at pc scale; and (ii) the extreme starburst at pc and kpc scale. This 
OF has generated 4 or more concentric expanding bubbles and the BAL systems. 

Specifically, inside and very close to the nucleus the 2D spectra show the presence 
of an OF emission bump in the blend Ha+[N n], with a peak at the same velocity 
of the main BAL-I system (VEj ec ti nBAL-i ~ -4700kms _1 ). This bump was more 
clearly detected in the area located at Of! 6-1'.' 5 (490-1220 pc), to the south-west of the 
nucleus core, showing a strong and broad peak. In addition, in the same direction (at 
PA ^60°, i.e. close to the PA of the small scale radio jet) at l'.'7-2'.'5, we also detected 
multiple narrow emission line components, with "greatly" enhanced [Nil] /Ha ratio 
(very similar to the spectra of jets bow shocks). These results suggest that the BAL-I 
system is generated in OF clouds associated with the parsec scale jet. 

The HST images show 4 (or 5) nuclear bubbles or shells with radius r ~ 2.9, 1.5, 
1.0, 0.6 and 0.2 kpc. For these bubbles, the 2D Ha velocity field (VF) map and 2D 
spectra show: (i) At the border of the more extended bubble (SI), a clear expansion 
of the shell with blueshifted velocities (with circular shape and at a radius r ~ 5'.'0). 
This bubble shows a rupture arc -to the south- suggesting that the bubble is in the 
blowout phase. The axis of this rupture or ejection (at PA ^00°) is coincident with the 
axis of the intermediate and large scale structures detected at radio wavelengths, (ii) 
In addition, in the 3 more external bubbles (SI, S2, S3), the 2D WHT spectra show 
multiple emission line components with OF velocities, of (VoFBubbic) SI, S2 and S3 
= [-(650-410) ±30], [-500±30],and [-230±30] kms" 1 . (iii) In all the circumnuclear 
region (1'.'8 < r < 5"), the [Nll]/Ha and [S n]/Ha narrow emission line ratios show 
high values (> 0.8), which are consistent with LINER/OF processes associated with 
fast velocity shocks. Therefore, we suggest that these giant bubbles are associated with 
the large scale nuclear OF component, which is generated -at least in part- by the 
extreme nuclear starburst: i.e., type II SN explosions. 

The variability of the short lived BAL-III Na ID system was studied, covering 
almost all the period in which this system appeared (between ~1984-2004). We found 
that the BAL-III light curve (LC) is clearly asymmetric with: a steep increase, a clear 
maximum and an exponential fall (similar to the shape of a SN LC). The origin of 
this BAL-III system is discussed, mainly in the frame work of an extreme explosive 
event, probably associated with giant SNe/hypernova explosions. 

Finally, the IR colours diagram and the UV-BAL systems of IR+GW/OF+Fe II 
QSOs are analysed. This study shows a new BAL IR QSO and suggest/confirm that 
© 2005 RAS these objects could be nearby young BAL QSOs, similar to those detected recently at z 

~ 6. 0. We propose that the phase of young QSO is associated with: accretion of large 
amount of gas (by the SMBH) + extreme starbursts + extreme composite OFs/BALs. 
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1 INTRODUCTION 

An important issue in astrophysics and cosmology is the 
study of extreme star formation and galactic winds/out 
flows, in mergers and QSOs and their relation to the early 
phases of the formation of galaxies and AGNs (see for ref- 
erences Lipari et al. 2004a,b,c,d). 



1.1 Luminous IR galaxies, mergers and AGNs 

Luminous and ultraluminous IR galaxies (LIRGs: Lib > 
10 11 L Q and ULIRGs: Li R > 10 12 L Q , respectively) are 
dusty, strong IR emitters where frequently a strong en- 
hancement of star formation is taking place (for references 
see Lipari et al. 2004a, b,c,d). Imaging surveys of LIRGs 
and ULIRGs show that a very high proportion ( ~ 70-95%) 
are mergers or interacting systems (Joseph & Wright 1985; 
Rieke et al. 1985; Sanders et al. 1988a; Melnick & Mirabel 
1990; Clements et al. 1996). In luminous IR galaxies there 
is a clear increase of the nuclear activity with the increase 
of the IR luminosity (Sanders et al. 1988a, b; Veilleux et al. 
1999, 2002). 

Strong evidence indicates that virtually all luminous 
IRAS galaxies are rich in interstellar gas, which is highly 
concentrated in their nuclei (Sanders et al. 1987, 1988a; 
Sanders, Scoville & Soifer 1991; Scoville et al. 1991). Specifi- 
cally, there are observational evidence and theoretical works 
suggesting that in IR mergers and IR QSOs tidal torque and 
loss of angular momentum drive large amount of interstel- 
lar gas into the central regions, leading to extreme starburst 
processes and probably fuelling a supermassive black hole 
(see for references Lipari et al. 2004a, d) 

An important, common phenomenon associated with 
the IR galaxies is the galactic wind: we found observa- 
tional evidence for galactic wind (GW) features from star- 
bursts and/or AGNs in luminous IR systems (Lipari et al. 
2004a,b,c,d, 2003, 2000a,b, 1994, 1993). 



1.2 Galactic winds (associated with starbursts 
and AGNs) 

Galactic winds and outflows have been observed mainly in 
starburst and Seyfert galaxies (see Heckman et al. 2000, 
1990; Cecil et al. 2002; Veilleux et al. 2002a, 1993; Lipari et 
al. 2004a,b,c,d, 2003, 2000b, 1994). There is substantial the- 
oretical literature about galactic winds associated with both 
processes: starbursts (Larson 1974; Ostriker & Cowie 1981; 
Chevalier & Clegg 1985; Ikeuchi & Ostriker 1986; Tomisaka 
& Ikeuchi 1988; Norman & Ikeuchi 1989; Suchkov et al. 1994; 
Strickland & Stevens 2000; and others) and with AGN (see 
for references Veilleux et al. 2002a). 

There are clear evidence of galactic -and local- winds, 
shells, arcs and bubbles generated by multiple SN explosion 
and massive star winds in starbursts. Our understanding of 
the main phases of galactic winds associated with starbursts 
was improved significantly by the use of theoretical and nu- 
merical models (see Strickland & Stevens 2000; Suchkov et 
al. 1994; Mac Low, McCray & Norman 1989; Tomisaka & 
Ikeuchi 1988). In general, good agreement has been found 
between these models and the observations. 

On the other hand, for galactic winds associated with 



AGN the situation is more complex, and very different mod- 
els are proposed in order to explain the observed data. In 
these models the OF could be generated by jets driven ther- 
mal winds, accretion discs winds, X-ray heated torus winds, 
etc. (see Veilleux et al. 2002a; Cecil et al. 2002; Capetti 2002; 
Morganti et al. 2003). 

However, there are increase evidence that nuclear galac- 
tic jets are one of the main sources of out flows in 
QSOs/AGNs (Cecil et al. 2002; Veilleux et al. 2002a; Capetti 
2002). In particular, there are clear detection and evidences 
of the interaction of jets with various phases of the ISM, 
including the narrow line region (NLR) . 

1.3 IR Mergers with galactic winds 

IR mergers often show strong starbursts with powerful out 
flow (OF) and galactic winds (GW; Lipari et al. 2003, 
2004a, b,c,d). Comparing our OF data base of IR merg- 
crs/QSOs (Lipari et al. 2004a, their Table 1) with two sam- 
ples of nearby IR mergers we found in both samples a high 
proportion of mergers with galactic winds: ~75%. These 
results suggest (or confirm) that: (i) GWs are "frequent 
events" in IR mergers, and (ii) extreme starbursts + GW 
and extreme IR emission could be simultaneous processes, 
induced by merger events. This last conclusion is also sup- 
ported by a clear trend found in the plot of OF velocity vs. 
log Lir (see their Fig. 13), in the sense that extreme OF ve- 
locities are detected only in extreme IR emitters (ULIRGs) . 
We note that the reported trend is observed mainly for IR 
mergers (Lipari et al. 2004a). An explanation for this ob- 
served trend is that high values of OF velocity and Lir are 
both associated mainly with the same process: "starburst 
+ QSO" events, probably induced by mergers (already sug- 
gested by Lipari et al. 2003; see also Genzel et al. 1998; 
Veilleux et al. 1999; Sanders et al. 1988a; Rieke et al. 1985). 

The low velocity OF process (LVOF, Vlvof < 
700kms _1 ; see Lipari et al. 2004a,b, 2003, 2000a,b) ob- 
served in NGC 3256, 2623, 4039, 5514, etc -using 2D 
spectroscopy- is consistent with those found previously in 
Arp 220, Mrk 266, NGC 1614, NGC 3690, and other IR merg- 
ers with massive starbursts. In addition, we found an inter- 
esting fact: starbursts and LINERs are the main sources of 
ionization in "low velocity OF" IR mergers (Lipari et al. 
2000b, 2003, 2004a). 

1.4 IR QSOs with galactic winds 

The discovery and study of IR QSOs (see for references 
Lipari et al. 2004d, 2003; Zheng et al. 2002) raises several in- 
teresting questions, in particular whether they are a special 
class of QSOs. We found, or confirm, that a high percent- 
age of IR QSOs show extreme out flows with giant galactic 
arcs, merger features, BAL systems, and extreme Fe II emis- 
sion, and are radio quiet (Lipari et al. 2004b, 2004a, 2003; 
Lipari, Terlevich & Macchetto 1993; Lipari, Macchetto & 
Golombek 1991a). We suggested that these objects could be 
young IR— active galaxies at the end phase of a strong star- 
burst: i.e., composite and transition QSOs (Lfpari 1994). 

Recently, for IR QSOs we found: (i) extreme velocity OF 
(EVOF, Vevof > 700 km s _1 ) objects with a composite nu- 
clear source: starbursts + AGNs/QSOs; for e.g. in Mrk 231, 
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IRAS 19254-7245, 01003-2230, 13218+0552, 11119+3257, 
14394+5332, and others (Li'pari et al. 2003); (ii) high res- 
olution HST WFPC2 images of IR + BAL + Fe n QSOs 
show in practically all of these objects arcs or shell features 
probably associated with OF or merger processes (Li'pari et 
al. 2003). 

Recently, using our data base of GW/OF in IR merg- 
ers/QSOs we found that in all the IR QSO candidates the 
H/3 broad line component is blushifted in relation to the 
narrow one, which is clearly consistent with the result ob- 
tained from the study of strong Fe II + BAL emitters, by 
Boroson & Meyer (1992). They proposed that blueshifted 
offset/asymmetry detected in the Ha broad components-of 
IR QSOs with strong Fe II + BAL systems-is probably due 
to the emission of the out flowing material, associated with 
the BAL process. Furthermore, our measurements are also 
consistent with the results of the study of multiple emis- 
sion line components, in some starburst nucleus of galaxies 
(Taniguchi 1987). Thus, a possible explanation for "at least 
part" of these blushifted broad line systems in IR QSOs is 
the high speed OF/GW generated in a burst of SN event 
near the nuclear region (Heckman et al. 1990; Terlevich 
et al. 1992; Perry & Dyson 1992; Lipari et al. 2003, 1994, 
2003, 2004a; Taniguchi et al. 1994; Scoville & Norman 1996; 
Lawrence et al. 1997; Collin & Joly 2000; and others). 

1.5 BALs in IR mergers/QSOs 

Some of the results obtained for nearby BAL QSOs, such 
as strong IR and Fen emission, strong blue asymmetry/OF 
in Ha, radio quietness, and very weak [O ni]A5007 emission 
(Low et al. 1989; Boroson & Meyers 1992; Li'pari et al. 1993, 
1994, 2003; Lipari 1994; Turnshek et al. 1997), can be ex- 
plained in the framework of the starburst+AGN scenario. 
In our study of Mrk231 and IRAS 0759+6559 (the nearest 
extreme IR + GW/OF + Fen + BAL systems), we detected 
typical characteristics of young-starburst QSOs. In our evo- 
lutive model for young and composite IR QSOs (see for ref- 
erences Li'pari 1994) suggested that BAL systems could be 
linked to violent supermassive starburst+AGN which can 
lead to a large-scale expanding shells, often obscured by 
dust. Several articles suggested that this evolutive model 
shows a good agreement with the observations (see Canal- 
izo & Stockton 1997; Lawrence et al. 1997; Canalizo et al. 
1998). 



Furthermore, the important detection of a correlation 
between the mass of galactic bulges and the mass of su- 
permasive black hole is a confirmation that the formation 
and evolution of galaxies (bulges/ellipticals, mergers) and 
supermassive black hole (AGNs and QSOs) are physically 
related to one another (Magorrian et al. 1998; Ferrarese & 
Merrit 2000; Gebhardt et al. 2000; Kormcndy 2000; Merrit 
& Ferrarese 2001; Kormendy & Richstone 1995). 

In the last years, several possible links between merg- 
ers, starbursts, IR QSOs and ellipticals have been proposed. 
Specifically, Joseph et al., Sanders et al. and Li'pari et al. sug- 
gested three complementary sequences and evolutive-links: 

(i) merger — > giant shocks — > super-starbursts + galac- 
tic winds — > elliptical galaxies; 

(ii) merger — > H2-inflow (starbursts) — > cold ULIRGs 
-> warm ULIRGs + QSOs; 

(iii) merger/s — > extreme starburst + galactic- wind (in- 
flow + outflow) — > IR + Fen + BAL composite/transition 
QSOs — * standard QSOs and ellipticals. 

Recently, several studies have confirmed the composite 
nuclear nature and the presence of merger and galactic wind 
features in several IR QSOs (Canalizo & Stockton 2001). In 
particular, Lipari et al. (2003, 2004a,d) found kinematical 
and morphological evidence of galactic winds in nearby "IR 
+ Fell + BAL QSOs". 

Thus, the study nearby "IR + Fell + BAL compos- 
ite/transition QSOs" is one observational way to analyse 
the possible relation between IR mergers and IR QSOs. We 
have started detailed morphological, spectroscopic and kine- 
matics studies of this kind of objects (since 1987: see Li'pari 
et al. 1991a, 1993, 1994, 2003, 2004d; Lipari 1994). 



1.6 Possible links between IR mergers and IR 
QSOs 

The luminosities and space densities of ULIRGs in the local 
Universe are similar to those of quasistellar objects (QSOs; 
Soifer, Houck & Neugebauer 1987). In addition, at the high- 
est IR luminosities, the presence of AGNs (and mergers) in 
LIRGs becomes important. Thus LIRGs probably represent 
an important stage in the formation of QSOs and elliptical 
galaxies. These results strongly suggest that it is important 
to perform detailed studies of possible links among mergers, 
ULIRGs, QSOs and elliptical galaxies (Toomre & Toomre 
1972; Larson 1974; Toomre 1977; Schweizer 1982; Joseph & 
Wright 1985; Sanders et al. 1988a,b; Li'pari et al. 1993, 1994, 
2003; Sanders & Mirabel 1996; Colina et al. 2001). 
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2 THE PROGRAMME AND OBSERVATIONS 
2.1 The programme 

2.1.1 IR mergers/ 'QSOs with galactic winds (the main 
programme) 

A current key issue in astrophysics is to explore the evolu- 
tion of the star formation process, especially at high redshift, 
when the galaxies/QSOs formed, and where it is expected 
that the star formation rate is very high. With these aims 
in mind we began a project to study nearby star forming 
+ GW galaxies and distant Lya emitters (see Lipari et al. 
2004a, b,c,d). The first step in this project is to understand 
the star formation process in nearby galaxies because we can 
obtain more detailed and unambiguous information. Thus, 
our groups started a study of nearby IR mergers/QSOs, 
which are an excellent laboratory-at low redshift-for the 
analysis of extreme star formation and GW processes (Li- 
pari et al. 2004a, 2003, 2000a,b, 1994). 

This study is based mainly on integral field spec- 
troscopy, obtained at the European Northern Observatory 
(ENO, La Palma-Spain) , the European Southern Observa- 
tory (ESO, Chile), the Complejo Astronomico El Leoncito 
(CASLEO, Argentina), and Bosque Alegre (BALEGRE, Ar- 
gentina) observatories, with the 4.2 m, 3.6 m, 2.15 m and 1.5 
m telescopes, respectively. The characteristics and goals of 
the programme have been described in detail by Lipari et al. 
(2004a). In particular, the main goals of this programme (at 
low redshift) is to analyse in detail the properties of the dif- 
ferent stages of extreme starbursts, galactic winds, mergers, 
QSOs, and elliptical galaxies (and their interrelation). 

In the present paper we present new results from our 
programme of study of "nearby" IR mergers/QSOs with 
galactic winds (EVOF): for Mrk231. 

2.1.2 BALs in IR + GW + Fell mergers/QSOs 

A study and search of UV (and optical) BAL systems in 
extreme IR + GW/OF + Fe n QSOs have been started 
since 1993 (see Lipari 1994; Zheng et al. in preparation), 
using mainly IUE and HST UV-spectra, plus ESO NTT, 
KPNO, APO and CASLEO optical data. In general, the 
role of BALs in IR+GW/OF+Fcll QSOs/mergers must be 
carefully considered, since: (i) Low et al. (1989) and Boroson 
& Meyers (1992) found that IR selected QSOs show a 27% 
low-ionization BAL QSO fraction compared with 1.4% for 
the optically selected high-redshift QSOs sample (Weymann 
et al. 1991); (ii) extreme IR galaxies (ULIRGs) are mainly 
mergers (Section 1); (iii) very recently Maiolino et al. (2003) 
reported also a high fraction of BAL QSOs at very high 
redshift (z ~6). The high percent of occurrence of broad 
absorption in extreme IR + GW/OF + Fell QSOs/mergers 
may be signals a fundamental relation (rather than merely 
a coincidence), and deserves detailed studies. 

Lipari et al. (1993, 1994, 2003, 2004a,d); Scoville & Nor- 
man (1996); Egami et al. (1996); Lawrence et al. (1997) and 
others proposed that the extreme IR + GW/OF + Fen 
-(- BAL phenomena are related -at least in part- to the 
end phase of an "extreme starburst+AGN" and the associ- 
ated "powerful bubble/galactic-wind". At the final stage of 
a strong starburst, i.e., type II SN phase ([8-60] x 10 6 yr from 
the initial burst; Terlevich et al. 1992; Norman & Ikeuchi 



1989; Suchkov et al. 1994) giant galactic arcs and extreme 
Fcll+BAL systems can appear. 

Specifically, the search of UV broad absorption line was 
focused mainly in the IR objects of our data base of IR + 
GW/OF + Fen QSOs/mergers (see Table 4), and studying 
the redshifted lines: Lya, Si IV, C IV, and Mg II. Lipari 
(1994) reported a first study of the UV BAL systems in 
IRAS 0759+651, with projected ejection velocity of about 
10000 kins -1 . Previously, it was also reported UV BAL sys- 
tems in Mrk 231 (IRAS 12540+5708), IRAS 17002+5153, 
and IRAS 14026+4341. 

It is important to note that in Table 4 several exam- 
ples clearly show the consistency of the OF data presented. 
In particular, for IRAS F05024-1941, our study of EVOF in 
the ULIRG 1 Jy Sample (using "emission" lines), we found 
that this object is a candidate to EVOF, with a velocity 
close to 1500 km s _1 (Lipari et al. 2000b). From a study 
of narrow "absorption" lines Rupke et al. (2002) obtained 
for this object an OF value of 1600 km s _1 . We have ver- 
ified that for Mrk 231 the results using the offset method 
(by Zheng et al. 2002) gave the same value of OF: -1000 
km s _1 , to that already obtained from the detection of two 
emission line systems in [O ll]A3727 (by Lipari et al. 1994). 
A similar agreement between the values of the OF obtained 
from multiple emission lines (Lipari et al. 2000b) and the 
offset method (Zheng et al. 2002) was found, for the cases 
of IRAS 11119+3257 and IRAS 15462-0450 (see Table 4). 

In this paper we report new results from this pro- 
gramme. In particular, the detection of a new IR + GW/OF 
+ Fe II QSO with BAL system. 

2.2 Mrk 231: the nearest extreme IR + GW + 

BAL + Fell merger/QSO (of the programme) 

Since Markarian 231 was discovered (by its strong ultraviolet 
continuum; Markarian 1979) it has been extensively studied 
and has been recognised as a remarkable galaxy: it is the 
most luminous galaxy in the local universe (z < 0.1), the 
nuclear spectrum shows extreme broad emission lines and 
strong BAL systems, the morphology is associated with a 
merger, etc (Arakelian et al. 1991; Adams 1972; Adams & 
Weedman 1972; Rieke & Low 1972, 1975; Boksemberg et al. 
1977; and others). In particular, it is important to remark 
that Sanders et al. (1987); Bryant & Scoville (1996) and 
Downes & Solomon (1998) found a very high concentration 
of CO molecular gas in the nuclear region (3 x 10 9 M Q , 
within a diameter of l'.'O, 814 pc). 

Before the IRAS was lunched, Mrk 231 was known to 
have an extreme IR luminosity with Mk = -24.7 (Rieke 
& Low 1972, 1975; Cutri, Rieke & Lebofsky 1984). Even 
after IRAS expand the known population of ULIRG, Mrk 
231 remain one of the most luminous object in the local 
universe, with L /fl [ 8 _ 1000Mm j = 3.56 x 10 12 L Q , and L IR /L B 
— 32. Consequently, the bolometric luminosity (dominated 
by the continuum IR emission) places the nucleus of this IR 
merger among the QSOs. Weedman (1973) already proposed 
that this extreme IR luminosity is associated with ultraviolet 
(UV) emission re-radiate thermally -in the IR- by dust. 

Mrk 231 also shows very interesting spectral character- 
istics, dominated in the optical by extremely strong Fe II and 
broad Balmer emission lines at Z em = 0.042 (Arakelian et 
al. 1971; Adam & Weedman 1972; Boksemberg et al. 1977). 
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In addition, it shows remarkable absorption line systems: a 
clear stellar absorption at Z a bs = 0.042 plus at least three 
strong broad absorption line (BAL) systems. These strong 
BAL systems show the following velocity of ejection: V e j ec t 
of BAL I, II and III -4700, ~6000, and ~8000 km s" 1 . The 
last BAL-III system sudden appear sometime between 1984 
December and 1988 May (Boroson et al. 1991) 

The extreme IR luminosity of Mrk 231 is associated 
with two main sources of energy: (i) an AGN detected clearly 
at radio and X-ray wavelengths (see for references Condon et 
al. 1991; Ulvestad et al. 1999; Turner 1999; Gallagher et al. 
2002), and (ii) an extreme nuclear + circumnuclear starburst 
observed almost at all wavelengths (see Downes & Solomon 
1998; Taylor et al. 1999; Gallagher et al. 2002). In addition, 
Lipari et al. (1994) reported kinematics, morphological and 
physical evidence of a powerful nuclear outflow process: i.e. 
a "galactic superwind with a broken shell, in the blowout 
phase". Even Lipari et al. (1994, 1993) suggested that Mrk 
231 is an evolving "young and composite IR QSO". 

Thus, Mrk 231 is the nearest extreme IR + GW/OF + 
BAL + Fen QSOs. In this paper, new results of the OF 
process in Mrk 231 are reported. Throughout the paper, a 
Hubble constant of Ho = 75 km s _1 Mpc -1 will be assumed. 
For Mrk 231 we adopted the distance of ~168 Mpc (cz = 
12465 ±10 km s _1 ), and thus, the angular scale is 1" ps814 
pc. 



2.3 La Palma/WHT + INTEGRAL 2D 

spectroscopy 

Two-Dimensional (2D) optical spectroscopy of Mrk 231 was 
obtained during a photometric night in April 2001 at the 
4.2 m telescope at Observatorio Roque de los Muchachos on 
the Canary island of La Palma. The 4.2 m William Herschel 
Telescope (WHT) was used with the INTEGRAL fibre sys- 
tem (Arribas et al. 1998) and the WYFFOS spectrograph 
(Bingham et al. 1994). The seeing was approximately l'.'O. 

INTEGRAL links the f/11 Naysmith focus of WHT 
with the slit of WYFFOS via three optical fibre bundles. A 
detailed technical description of INTEGRAL is provide by 
Arribas et al. (1998); here we only recall its main character- 
istics. The three bundles have different spatial configuration 
on the focal plane and can be interchanged depending the 
scientific programme or the seeing conditions. At the focal 
plane the fibres of each bundles are arranged in two groups, 
one forming a rectangle and the other an outer ring (for 
collecting background light, in the case of small-sized ob- 
jects). The data analysed in this paper were obtained with 
the standard bundle 2: which consistent of 219 fibres, each 
of 0'.'9 in diameter in the sky. The central rectangle is formed 
by 189 fibres, covering a field of view of 16'.'4 x 12'.'3; and 
the other 30 fibres form a ring 90" in diameter. 

The WYFOS spectrograph was equipped with a 1200 
lines mm -1 grating, covering the AA6200-7600 Arange. A 
Tex CCD array of 1124 x 1124 pixels of 24 /im size was used, 
giving a linear dispersion of about ~1.45 A pixel -1 (~2.5 A 
effective resolution, ~ 70 km s _1 ). With this configuration 
and pointing to the central region of Mrk 231, we took three 
exposure of 1800 s. 



2.4 HST ACS, WFPC2 and NICMOS broad band 
images 

Near-UV HST Advanced Camera for Surveys (ACS) 
archival images of Mrk 231 were analysed, obtained with 
the High Resolution Channel (HRC). They include images 
with the filter F330W (3354 A, AA 588 A). The scale is 
0"027 pixel -1 . 

Optical HST Wide Field Planetary Camera 2 (WFPC2) 
observations were analysed, which include broad-band im- 
ages positioned on the Planetary Camera (PC) chip with 
scale of 0'.'046pixer 1 , using the filters F439W (4283 A, AA 
464 A, ~B Cousin filter), and F814W (8203 A, AA 1758 A, 
-I). 

HST Near IR Camera and Multi Object Spectrometer 
(NICMOS) archival data were studied, which include mainly 
images with the filter F160W (1.60 /im, AA 0.4 /im), using 
the camera 2 with a scale of 0'.'076 pixel -1 . 

For details of the HST observations see also Table 1. 

2.5 HST UV-FOS and UV-STIS ID Spectroscopy 

HST FOS aperture spectroscopy of the main optical nu- 
cleus of Mrk 231 was obtained, from the HST archive 
(at ESO Garching). The spectra were taken with the 
G190L (AA1150-2300), G190H (AA1575-2320) and G270H 
(AA2225-3290) gratings and the blue detector. The G190H 
and G270H observations were made with the effective aper- 
ture of 4'.'3 x l'/4; and these spectra have a resolution of 
~3 and 4 A, FWHM respectively. The G160L observations 
were made with the aperture of l'.'O, and the spectra have a 
resolution of ~8, FWHM. 

In addition, HST FOS spectra of extreme IR + Fe II 
QSOs with GW were analysed. These data were used in our 
study/search of UV-BAL systems (for details see Table 1). 

HST STIS long slit spectra of IR + GW/OF + Fe II 
QSOs were obtained (mainly from the HST archive). The 
spectra were taken with the G750M grating, with the slit 
07 1 x 5270, giving a resolution of 50 km s -1 . The STIS 
CCD has a scale of 0'.'05 pixel" 1 . 

2.6 La Palma/NOT optical broad band images 

Optical V broad band optical image of Mrk 231 was ob- 
tained during a photometric night in 1991 May at the 2.5 m 
Nordic Optical Telescope (NOT) at Observatorio Roque de 
los Muchachos on the Canary island of La Palma. A TEX 
CCD chip with scale 07197 per pixel was used. The obser- 
vations were carried out in excellent seeing conditions, with 
a values of 077 (FWHM). 

2.7 KPNO long-slit ID spectroscopy 

Optical long slit spectra of Mrk 231 and IR + GW/OF + 
Fe II QSOs were obtained on the KPNO 2.15 m telescope 
with the Gold Cam Spectrograph, by using two grating 
of 300 grooves mm -1 (7 A ~300 km s _1 resolution) cover- 
ing the wavelength ranges 3350-5200 and 5100-7100 A. The 
observations were taken during one photometric night in 
February 1991 (Table 1). The spectra were obtained with 
a slit width of 1.5 arcsec. The seeing was in the range 1.0- 
1.5 arcsec (FWHM). 
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The data of the IRAS 1 Jy ULIRGs sample were ob- 
tained with the Gold Cam Spectrograph on the KPNO 
2.15 m telescope by using a grating of 300 grooves mm _1 
(8 A resolution); and covering the wavelength range 4500- 
9000 A. 

2.8 IUE UV ID spectroscopy 

UV International Ultraviolet Explorer (IUE) spectra were 
obtained for our programme of search of UV-BAL systems 
in extreme IR + GW/OF + Fe n QSOs (Table 1). In ad- 
dition, data from the IUE archival -of this type of objects- 
were analysed. The IUE LWP and SWP spectra cover the 
wavelength ranges 1150-1980 and 1900-3290 A,; with a res- 
olution of ~8 A (~300 km s _1 ). 

2.9 ESO NTT and 3.6 m telescope spectra 

The ESO Faint Object Spectrograph and Camera II 
(EFOSC II) on the 3.5 m new technology telescope (NTT) 
at La Silla was used, to obtain long slit spectra of IR + 
GW/OF + Fe ii QSOs. These long-slit observations were 
taken during 1 photometric night in July 1993 (see Table 
1). The spectra cover the wavelength ranges ~3560-5200 A, 
with a resolution of ~7 A (~300 km s _1 ) 

The ESO Faint Object Spectrograph and Camera 
(EFOSC) on the 3.6 m telescope at La Silla was used to ob- 
tain long-slit spectra and high resolution images. Medium- 
resolution spectra were obtained with the B150, O150, and 
R150 grisms, which provide a resolution of ~ 7 A in the 
wavelength range 3600-9800 A (during 3 photometric nights 
in 1995 January and 1991 July; see Table 1). 

2.10 Apache Point Observatory (APO) ID 
spectroscopy 

The APO spectra of IR GW/OF + Fe II QSOs were taken 
with the 3.5 m telescope. Optical long-slit observations were 
obtained during 1 photometric night in September 2004 
(see Table 1). The medium resolution spectra give an effec- 
tive resolution of ~6 A (~290 km s _1 ), covering the wave- 
length ranges AA3500-5600 and 5300-10000 A. The seeing 
was ~1'.'5 (FWHM). 



an effective resolution of ~8 A (~300 km s x ) covering the 
wavelength range AA4700-7200. The seeing was in the range 
l'.'5-2'.'5 (FWHM). 

2.12 Reduction 

The IRAF* software package was used to reduce and anal- 
yse the data. The reduction of the 2D spectroscopic observa- 
tions consist of two main steps: (1) reduction of the spectra 
(for each 219 fibres), and (2) generation of 2D maps by inter- 
polated the selected parameter (e.g., emission- line intensity, 
continuum intensity, radial velocity, etc) from the grid val- 
ues defined by the fibre bundle. The step (1) was basically 
done in the same way as for long-slit spectroscopy, includ- 
ing bias subtraction, aperture definition and trace, stray- 
light subtraction, the extraction of the spectra, wavelength 
calibration, throughput correction and cosmic-ray rejection. 
We obtained typical wavelength calibration errors of 0.1 A, 
which give velocity uncertainties of 5 km s_i. For step (2) 
we used the software package Integral^ , with 2D interpola- 
tion routines. In particular we transformed ASCII files with 
the positions of the fibres and the corresponding spectral 
features, into regularly spaced standard FITS files. Maps 
generated in this way are presented in the next sections. 

The IRAF and STSDAS* software packages were used 
to analyse the HST, ESO, IUE, ESO, KPNO, NOT, APO, 
CASLEO imaging and spectrophotometric data. 

The HST data were reduced at the Space Telescope Sci- 
ence Institute (STScI), using the standard process. All the 
bands were calibrated according to the procedures detailed 
in Lipari et al. (2004a). Specifically, photometric calibration 
of the HST data was performed using published photometric 
solution (Lipari et al. 2004a). 

For the study of the kinematics data we used the soft- 
ware ADHOC§ 



2.11 CASLEO ID spectroscopy 

Spectrophotometric observations of IR + GW/OF + Fe II 
QSOs were taken at CASLEO (San Juan, Argentina) with 
the 2.15 m Ritchey-Chretien telescope. Optical long-slit and 
aperture spectroscopy observations were obtained during 3 
photometric nights in June 1989, July 1993, March 1997 (see 
Table 1). Long-slit spectroscopic observations with medium 
resolution were obtained with the University of Columbia 
spectrograph (UCS; e.g. Lipari et al. 1997). The medium 
resolution spectra were made using a 600 line mm -1 grating, 
a slit width of 2'.'5, which gives an effective resolution of ~6 A 
(~290 km s _1 ) and a dispersion of 120 A mm -1 , covering the 
wavelength range AA4000-7500 A. Aperture spectroscopic 
data were obtained with the 'Z-machine' (e.g. Lipari et al. 
1991a, b). These aperture spectra were made using a 600 
line mm -1 grating, giving a dispersion of 130 A mm -1 and 



* IRAF is the imaging analysis software facility developed by 
NOAO 

t Integral is the imaging and spectroscopic analysis software fa- 
cility developed by the Instituto de Astrofisica de Canarias (IAC) 
t STSDAS is the reduction and analysis software facility devel- 
oped by STScI 

§ ADHOC is the analysis 2D/3D kinematics software facility de- 
veloped by Marseille Observatory 
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3 RESULTS 

This section focuses on presenting (for Mrk 231): (i) high 
resolution HST near-UV, optical and near-IR images of the 
main body; (ii) 2D spectroscopy mainly for the central re- 
gions; (iii) ID spectroscopy for the nuclear region; (iv) a 
study of variability of the Na ID BAL-III system (based in 
ID spectroscopy). In particular, we study the extreme nu- 
clear OF with multiple expanding concentric shells, giant 
shocks and BAL systems. In addition, new observations of 
UV-BAL systems in IR+GW/OF+Fc n QSOs are analysed. 

3.1 Near-UV, optical and near-IR HST and NOT 
images: multiple concentric nuclear bubbles 
(and knots) 

Here, a multiwavelength morphological study (using mainly 
high resolution HST near-UV, optical and near-IR broad 
band images) of the nuclear and circumnuclear structures of 
Mrk 231 is presented. 

The NOT V image (Fig. 1) shows almost the whole 
merger. This galaxy consist of a nearly elliptical main body 
(R ~ 10 kpc) with a compact nucleus, plus two faint tidal 
tails (see also Lipari et al. 1994; Hamilton & Keel 1987; Neff 
& Ulvestad 1988). These are typical features observed in 
advanced mergers. 

It is important to note, that previously, in the circum- 
nuclear region of Mrk 231 Lipari et al. (1994) found clear 
evidence of a powerful nuclear galactic-wind and they pro- 
posed that this OF process (in part associated with a star- 
burst) has generated the circumnuclear blue arc/shell (at 
r ~ 3 kpc, and to the south of the nucleus). However, Ar- 
mus et al. (1994) suggested that this arc is associated with 
an obscured second nucleus (they also suggested that in this 
blue region or "shell" there is no evidence of star-formation 
process). HST-WFPC2 observations of Mrk 231 confirmed 
that this blue arc is a "dense shell of star-forming knots" 
(see Surace et al. 1998; Lipari et al. 2003). In the next sec- 
tions we will expand these previous studies. 

3.1.1 The multiple concentric nuclear bubbles / shells 

Figs. 2, 3 and 4, show four (or five) concentric supergiant 
galactic shells or bubbles with centre in the nucleus and 
showing several bright knots. These figures present high res- 
olution HST broad-band images obtained in the near-UV, 
optical and near-IR wavelengths through the filters F330W 
(~ U), F439W (~ B), F814W (~ /), and F160W (~ H). 
We note that Figs. 3a,b,c,d were obtained from the subtrac- 
tion of a smooth image (for each filter) of the main body of 
Mrk 231, in order to study the morphological residuals. 
These figures show intereting structures, in particular: 

(i) Bubble or Shell-Sl: This is the more external and well 
defined shell detected in Mrk 231, with a radius Rsi = 
3'.' 5 = 2.9 kpc. Lipari et al. (1994) already noted that this 
shell shows a clear elongation to the south (of the nucleus). 
At the border of this elongation there is a blue arc with sev- 
eral knots. They associated this arc with typical structures 
generated in the blowout phase of a nuclear galactic wind 
(with giant shocks) . These shocks probably generated a new 
extranuclear star formation process, with knotty structure. 



Fig. 2 shows interesting details of this bubble + arc in the 
HST ACS, WFPC2 and NICMOS images. In particular, the 
HST WFPC2/B image depicts clearly the complete shell SI. 
However the HST ACS/U image shows remarkable details of 
the blue compact knots located at the border of the southern 
arc. 

(ii) Bubble-S2: The next internal bubble, at R S2 = l'.'8 = 
1.5 kpc, shows also a well defined circular shape. However, 
this more internal bubble is clearly obscured in the south- 
east region by reddening, associated with the presence of 
dust. Again, the HST WFPC2/B image depicts more clearly 
the partial shell S2. It is important to remark that this shell 
shows a very interesting knot-2 in the south-west direction 
(see Table 2 and the next sections for details about the prop- 
erties of this knot). 

(iii) Shell-S3: The next bubble, with a radius Rs3 = 
l'/2 = 1.0 kpc, depicts also a circular shape, and their emis- 
sion is strong mainly in the south area. This structure is 
well defined in the HST ACS/U image, which shows also a 
knotty structure. This shell depicts a bright knot-3, to the 
south-east. 

(iv) Shell-S4: This shell is clearly observed mainly in the 
HST ACS/U image (Fig. 2a; we note that this ACS im- 
age has the best spatial resolution) and also in the ACS/U, 
WFPC2/B residual images (Figs. 3a,b). The radius of this 
shell S4 is R S i = 077 = 0.6 kpc. 

(v) Possible shell-S5: This possible shell is only observed 
in the HST optical (I-B) colour image (Fig. 4) and very close 
to the nucleus, with a radius Rss = 0'.'3 = 0.2 kpc. 

(vi) Close to the nucleus, the HST images show radial 
filaments. These structures are more clear in the near-IR 
image since at these wavelengths the dust absorption is sig- 
nificantly smaller. In particular, it was found that these fila- 
ments are more larger and clear in the direction north-south 
(PA ~00°) and at the position angle of the radio jet (PA 
~65°; see sections 4.2 and 4.3 for details about this jet). We 
note that these PAs are coincident with the directions of the 
two main OF processes in Mrk 231 (see the next sections). 
In our programme of study of galactic winds in IR merg- 
ers/QSOs, similar filaments were detected associated with 
the OF processes of NGC 3256, NGC 2623 and NGC 5514. 

This type of circular and concentric "bubbles or shells" 
could be associated mainly with symmetric explosive events. 
In sections 3.2 and 4.3 the kinematics, age, energy and pos- 
sible origin of these structures will be discussed. 

3.1.2 The circumnuclear knots 

Table 2 presents the location and properties of the main 
knots, in the nuclear and circumnuclear areas of Mrk 231. 
In addition, the positions of these knots are labeled in Fig. 
3b. Surace et al. (1998) already studied in detail these knots 
using HST photometry, in the B and I bands. In this pa- 
per, we will study mainly their properties in relation to the 
process of OFs and their location in the shells. 

From the data of Table 2, it is clear that in the south- 
ern part of the shell SI, the knots 5, 6, 11 and 12 show blue 
colours, probably associated with new star formation events. 
In particular, these knots emit a strong UV continuum, due 
to the presence of a large number of massive O and B stars. 
These massive stars are usually responsible for the rupture 
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phase of the galactic bubble and also are the progenitors of 
core-collapse SNe. In the supergiant bubble of NGC 5514 
we also detected the presence of several knots with massive 
WR stars (inside the areas of rupture of the external shell) . 
The radius of this bubble in NGC 5514 -just detected in the 
rupture phase- is also of ~3 kpc, i.e. the same value of the 
shell SI in Mrk 231. We already noted that the results of hy- 
drodynamical models for galactic winds (associated mainly 
with extreme starbursts) suggest that 3 kpc is the typical 
value of radius for the beginning of the blowout phase (after 
^8 Myr, from the initial starburst; Suchkov et al. 1994). 

The remaining knots show different colours (see Table 
2). Some of the knots depict red colours, which could be ex- 
plained -in part- by the presence of dust in star formation 
areas. In particular, the HST NICMOS image shows that 
the knot 14 is bright at the H-band, suggesting that the red 
colour in this knot (located in the southern arc of the shell 
SI) is due mainly to reddening by dust. However, for another 
two knots with very red colours, the 2D spectroscopy suggest 
another origin for the red colours. For the knot 2 (located in 
the shell S2) the spectrum shows properties typical of bow 
shocks, similar to those found in radio jets (see for details 
the next section). In addition, the knot 29 (positioned at the 
north-west area of the shell SI) also shows very red colour, 
and it is located in another OF region. Where the [S ll]/Ha 
and [N ii]/Hq maps depict extended radial filament, proba- 
bly emerging from the bubble (see section 3.4). For this area 
also Krabbe et al. (1997) reported OF ejection (of ~1400 km 
s _1 in Paa; from a study of 2D near-IR spectroscopy). 

In the next section, the kinematics and physical prop- 
erties of the multiple concentric shells and the circumnu- 
clear knots will be analysed (using mainly WHT 2D spec- 
troscopy) . 

3.2 Mapping the outflows in the nucleus and in 
the multiple supergiant bubbles (2D Spectra) 

The Ha+[N ii] emission lines profiles obtained simultane- 
ously with INTEGRAL for the central region of Mrk 231 
(16"xl2"; 13.0 kpcxlO.O kpc) are shown in Fig. 5. These 
are reduced data, but represent only a small part (~100 A) 
of the full spectral coverage (6200-7600 A) at each of the 189 
fibre location. The nucleus was posioned close to the centre 
of the field (at fibre 107), where it is evident the presence of 
broad emission lines. 

First, the OF at small scale will be analysed, in the 
region close to the nucleus (R < 1-2 kpc). In particular, it 
is important to remark the following main results: 

(i) In the fibre 101 (adjacent to the fibre of the nucleus) 
the 2D spectrum shows the presence of an strong OF emis- 
sion bump in the blend Ha+[N n], with the peak corre- 
sponding to the same velocity of the main BAL-I system 
(VEjoctionBAL-i ~ -4700 kms -1 ). 

This bump was detected in the area covered by this fibre, 
at 0'.'6-l'.'5 (490-1220 pc), to the south-west of the nucleus 
core at PA ~60° , showing a strong and broad peak (see Fig. 
6b). 

It is important to note that the direction/PA of fibre 101 
(from the main nucleus) is the same that the PA of the small 
scale radio jet (PA je t -60-70°). 

(ii) Following the same direction -at PA ~60°- in the 



next fibre (91; Fig. 6c) we found multiple narrow emission 
line components, with "greatly" enhanced [Nil] /Ha ratio, 
very similar to the spectra of jets bow shocks found in the 
Circinus galaxy (by Veilleux & Bland-Hawthorn 1997). This 
feature was detected at 1'.' 7-27 5 (1340-2030 pc), inside of the 
shell S2 and coincident with the knot 2 

(iii) For the nucleus (i.e., fibre 107, r <0.45"~ 366 pc) the 
blue OF bump was observed, but with very weak emission 
(Fig. 6a). 

Therefore, these results are consistent with the interest- 
ing fact that the blue OF Ha bump plus the BAL-I system 
could be related with OF clouds associated with the small 
scale jet. A discussion of this point will be presented in sec- 
tion 4.3. 

Now, the OF kinematics of the multiple shells will be 
analysed, for the central region (R < 6 kpc). In the 3 ex- 
ternal bubbles SI, S2 and S3 the 2D WHT spectra show 
multiple emission line components (see Fig. 6), specially in 
the bright knots. Where the following OF velocities values 
were measured: 

(i) For the shell SI, 

• to the south, in the knot 14 (fibre 065), 
(VoFShcii-si/South-Foes) = [-420 ±30] km s -1 

• to the east, in the knot 4 (fibre 084), 
(VbFShdi-si/Ea S t-F084> = [-485 ± 30] kms -1 

• to the north-west, in the knot 29 (fibre 131), 
(VoFShcii-si/NorthWcst-Fm) = [—500 ± 30] kms 1 

• to the west, in the fibre 077, 

(VoFShcii-si/Wcst-F077> = [-650 ± 30] kms -1 

(ii) For the shell S2, in the knot 2 (fibre 091), 
(VoFShoii-S2/F09i) = [-500 ± 30] kms -1 , and 

(iii) For the shell S3, in the knot 3 (fibres 103, 102) 

<V O FSheii-S3/Fi03> = [-230 ± 40] kms" 1 . 

These values are consistent with expanding bubbles, 
and their properties and origin will be discussed in section 
4.3 

3.3 Velocity field of the ionized gas, in the nuclear 
and central regions (2D Spectra) 

In order to study the kinematics of the ionized gas, in the 
central region of Mrk 231, we measured the velocities from 
the centroids of the stronger emission lines Ha, [N n]A6584 
and [S ii] AA6717-6731, fitting gaussians (with the software 
SPLOT and SPECFIT; for details see section 2). In Mrk 231, 
first we study the main emission line component; which is 
strong in all the central region. In the shells the presence of 
multiple OF components required a detailed study. 

Fig. 7 shows -for the ionized gas- the Ha velocity field 
for the central region ~16"xl2" ~13 kpcxlOfcpc, with 
a spatial resolution/sampling of 079. This map was con- 
structed using the techniques described in section 2, and 
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for the main component of the Ha emission lines. The er- 
rors vary from approximately ±10 km s _1 in the nuclear and 
central regions (where the emission lines are strong), to ~20 
km s _1 for the weakest lines away from the nuclear areas. 
Using the 2D Ha velocity field (VF) we have obtained a 
mean value of the systemic velocity Vsyst = 12645±10 km 
s _1 , this value was defined as the zero of the VF. 

The Ha isovelocity colour map (Fig. 7) shows interest- 
ing structures: 

(i) Coincident with almost the border of the more ex- 
tended bubble or shell (SI), parts of a ring with blueshifted 
values of velocities and circular shape was detected. We note, 
that this structure reach the south border of the field, and 
probably part of this ring (to the south) is located outside 
of our INTEGRAL field. The radius of this structure/ring 
is r ~ 5'.'0 ~ 4 kpc. 

(ii) Close to the centre of this ring (which it is located 
near to the nucleus) there is a small and circular lobe with 
blueshifted velocities. 

In addition, inside of the ring 4 symmetric redshifted lobes 
were also detected. 

(iii) To the north-west, part of an elongated blueshifted 
lobe was found. In this structure we have measured high 
values of velocities of A V ~-200 km s _1 . 

These results are clearly consistent with those found 
in section 3.2. In particular, the detection of a blueshifted 
ring coincident with the position of the shell SI confirm the 
OF/expansion of this shell (even in the main component of 
the emission lines). The presence of an elongate blueshifted 
structure to the north-east (of the INTEGRAL VF) is also 
consistent with the detection -in this area- of filaments in 
the [S ll]/Ha and [N n]/Ha maps (see section 3.4) and the 
previous detection of OF ejections (Krabbe et al. 1997). 



3.4 Mapping the emission line ratios and the 
ionization structure: large scale galactic 
shocks 

The set of 2D spectroscopic data, which cover the main 
structures of the central region of Mrk 231, allows the inves- 
tigation of the ionization structure and the physical condi- 
tions in the gaseous phases. Figs. 8 (a) and (b) show 2D maps 
(of 16'.'4 x 12'.'3, with 0'.'9 spatial sampling) of the emission 
line ratios [N n] A6583/Ha, and [S n] A6717 + 31/Ha. These 
maps were constructed using the techniques described in sec- 
tion 2 and are based on the main component of the emission 
lines. Figs. 8(c) and (d) show the values of these ratios for 
each fibre (which were used in the generation of the maps). 

Figs. 8(a) and (b) show interesting features. We note 
the following: 

(i) For the area where we measured the narrow emis- 
sion line components, i.e. in almost all the circumnuclear 
region (l'/8 < r < 5"; i.e. including the 2 more external 
shells/bubbles SI and S2), the [Nn]/Ha emission line ratios 
show high values (> 0.8). This result, obtained from the 2D 
map (Fig. 8a) was verified using the individual spectrum of 
each fibre (Fig. 8c). 

(ii) The [S n] A6717 + 31/Ha map (Fig. 8b) also shows 
high values in the circumnuclear areas. In particular, 



we found a range of values for this ratio: 3 > [S 
ii] A6717 + 31/Ha > 0.8 (see Figs. 8d,b). 

(iii) Specially, to the north-west the [S ii] A6717 + 31/Ha 
and [N ll]/Ha maps show filaments, which are probably as- 
sociated with ejection from the bubble. Lipari et al. (2004d) 
already discussed in detail the fact that the [S ll]/Ha map 
is one of the best tracer of filaments associated with OF and 
shock processes. 

Thus, in almost all the circumnuclear and central re- 
gions of Mrk 231 (including the shells) these INTEGRAL 
emission line ratio maps show high values (> 0.8), which 
are consistent with an ionization process produced mainly 
by shock-heating in a in the outflowing gas of the expand- 
ing shells (Lipari et al. 2004a, d; Dopita & Sutherland 1995; 
Heckman 1980, 1996; Heckman et al. 1987, 1990; Dopita 
1994). Furthermore, these ratios (for the main emission line 
components) are located in the area of fast shock velocities of 
~400-500kms _1 in the upper part of the [N n] A6583/Ha 
vs. [S ll]A6717 + 31/Ha diagram (published by Dopita & 
Sutherland 1995, their Fig. 3a). 

Similar results were obtained in the 2D studies of the 
bubble NGC 3079 (Veilleux et al. 1994) and the OF neb- 
ula and bubble of NGC 2623 and NGC 5514 (Lipari et al. 
2004a,d). They found that the ratios [N ll]A6583/Ha and [S 
ll]A 6717 + 31/Ha are > 1, in almost all the bubbles and the 
OF regions. They associated these results to the presence of 
large scale OF + shocks. 

3.5 The optical BAL systems and the nuclear 
spectrum of Mrk 231 (ID— spectra) 

3.5.1 Optical BAL systems: the variability of the Na ID 
BAL III systems 

It is important to study in detail the optical absorption and 
emission features of this IR QSOs, in order to detect similar- 
ities and differences between this object and similar extreme 
IR + GW/OF + Fc II emitters. The optical and UV spec- 
trum of Mrk 231 (see Fig. 9a) presents the features that 
are typical of extremely strong optical Fe II emitters with 
BALs: i.e., strong Na I DAA5889-5895 absorption, strong 
unresolved blends of optical Fe II lines, and very weak high- 
excitation forbidden emission lines. It is important to remark 
the strong fall in the continuum flux at the UV wavelengths 
(this is also a typical feature of dusty luminous infrared 
galaxies; see section 3.6). 

Boksemberg et al. (1977) detected two nonstellar ab- 
sorption lines systems (I and II) with velocities of 6250 and 
8000 km s _1 , respectively. This correspond to ejection veloc- 
ities of -6650 and -4950 km s _1 , with respect to the systemic 
velocity of 12900 km s" 1 . Boi 'oson et al. (1991) detected a 
new absorption system III in the Na ID Ca H,K and He I, in 
their optical spectra of 1988. This system was not present in 
the spectra of Boksemberg et al. (1977), taken in 1975. With 
an observed radial velocity of 4660 km s~ , this new sys- 
tem has the highest ejection velocity -8240 km s _1 . The ori- 
gin of these absorption-line systems have been discussed by 
Boroson et al. (1991); Kollatschny, Dietrich & Hagen (1992); 
Boroson & Meyers (1992); Lipari et al. (1994); Foster, Rich 
& McCarty (1995); Smith et al. (1995); Rupke et al. (2002), 
and others. 

It is interesting to recall about the detection of another 
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absorption system (IV), which correspond to the systemic 
velocity and has the narrower absorption width. This system 
was already associated by Boksemberg et al. (1977) with hot 
stars in the nucleus of Mrk 231. 

Fig. 9b shows the BAL systems at the line Na ID. 
The corresponding width of these absorptions systems are: 
FWHM of BAL system I, II and III = 700, 130, 80 km s _1 , 
respectively. These are in general low values for BAL sys- 
tems. Thus, these values are more consistent with those of 
mini-BALs or associated absorption lines (AAL; de Kool et 
al. 2001, 2002). 

Now, we will expand the previous study of variability 
of the Na ID BAL III system (covering almost all the period 
in which this system appeared). In order to study the be- 
haviour of the variability of this absorption system, the mea- 
sured values of the equivalent width ratio of BAL III/(BAL 
I + II) were included in Table 3. These data were obtained 
from our observation and also derived from published spec- 
tra (together with data previously published by Kollatschny 
et al. 1992; and Forster et al. 1995). We note that the pre- 
vious study of variability cover the period 1980-1992 (with 
only few point after the maximum) . The observation of Ta- 
ble 3 cover the period between ~1980-2000. Figs. 10a,b show 
the shape of the BAL-III light curve (LC), which is clearly 
asymmetric with: a steep increase, a clear maximum and an 
exponential fall. In particular Fig. 10b shows a very good 
fit of the LC fall, using an exponential function. In general, 
the shape of the BAL III system LC is similar to those LC 
of SNe flux, emission line, etc. 

The probable origin of this BAL III system will be dis- 
cussed in section 4.3. Specifically, we will discuss 4 possible 
scenarios. One of this scenario is an explosive event asso- 
ciated with explosions of giant SNe/hypernova (with very 
massive progenitors and located close to the AGN and/or in 
accretion disks of AGNs). 

It is important to note that previously it was assumed 
that this light curve (of BAL III) is symmetric. This type 
of variability was associated mainly with absorbing clouds 
crossing the central continuum light source (Kollatschny et 
al. 1992). 



3.5.2 Out flow from the nuclear optical emission line 

For the nuclear region of Mrk 231, a new detailed study of 
multiple emission line components was performed. In par- 
ticular, for the line [O ll]A3727 the presence of at least 3 
components were detected (see Fig. 9c): with OF1 and OF2 
velocities of -800 and -1300 km s" 1 . The value of the OF 
previously reported by Lipari et al. (1994) was obtained from 
the blend of the OF1+OF2. In addition, in the line Ni II an 
OF component was measured, and it was measured a value 
similar to the OF1, of the line [O n]. 

It is interesting to note that Schmidt & Miller (1985) 
from a spectropolarimetry study detected an OF velocity 
(for the dust) of 700 km s~ . Which is a value very close 
-within the errors- to our OF1 velocity. 

A discussion of the extreme optical Fe II emission of 
this IR QSO, was reported by Boksemberg et al. (1978) 
and Lipari et al. (1994). A ratio of Fe Hopt/U/3 ~8, [O 
III5007/H/3 ~0.01, Na ID/Ha ~ 0.03 and Ha/H/3 ~5 were 
measured. These results are in agreement with the typical 



values obtained for extreme Fe IIopt emitters, like IRAS 
07598+6508, PC 1700+518, IRAS 18508-7815 (Lipari 1994). 



3.6 The UV BAL systems and spectrum of Mrk 
231 and IR + GW/OF + Fe II QSOs 

3.6.1 UV spectrum and BALs in Mrk 231 

Significant reddening in luminous infrared sources hampers 
the study of their UV properties. For Mrk 231, Figs. lla,b 
and 9a show -at UV wavelengths- a strong reddening, with 
weak Mg II, C IV and Lya emission lines. UV absorption as- 
sociated with the BAL system I are present in the Mg II and 
C IV lines (Figs. lla,b). In particular, The low-ionization 
emission line Mg II is observed superposed to the BAL sys- 
tem I (Fig. 11a). We measured for this absorption line a 
width (FWHM) of - 700 km s" 1 . This value is equal to the 
that measured -in this system- in the Na ID optical line. In 
addition, in the region between AA2000-2800 there are weak 
absorption due to UV Fe II. 

It is important to study in detail the UV (plus opti- 
cal) continuum of extreme IR-Fe II emitters and specifically 
BAL-QSOs because they show a clear deficit in the UV flux 
as compared to the non-BAL QSOs (see Weymann et al. 
1991; Sprayberry & Foltz 1992). Fig. 9a shows a clear drop 
in the level of the UV (and even in the B) continuum of Mrk 
231, for A shorter than 4500 A. We found that an SMC-like 
reddening of E(B-V) ~ 0.13 accounts for the differences be- 
tween the UV spectrum of Mrk 231 and the non-BAL QSO 
composite/mean spectrum (Reichard et al. 2003). 



3.6.2 UV-BALs in IR + GW/OF + Fe n QSOs 

From our programme of study and search of UV and optical 
BAL systems in extreme IR + GW/OF + Fell QSOs we 
already found interesting results in relation with the UV 
low ionization BAL system detected in IRAS 07598+6508 
(Lipari 1994). 

We recall that the standard definition of BAL QSO 
(Weymann et al. 1991) is based in the measurement of the 
equivalent width of the C IV resonance absorption line sys- 
tem (called balcity index: BI). The QSOs with BI > km 
s _1 are considered BAL QSO. In part by this reason we ob- 
served the BAL candidates at UV (for C IV and Mg II lines) 
and optical (for Na ID absorption) wavelengths. 

From this programme, the following main results were 
found: 

(i) The XJV-IUE spectrum of IRAS 21219-1757 (Figs. 
12a, b) clearly shows a blend of emission and absorption fea- 
tures -at low S/N- in the C IV line. In these two figures 
the UV spectrum of IRAS 21219-1757 is superposed with 
the spectra of standard low ionization BAL and non-BALs 
QSO, respectively (obtained from Reichard et al. 2003). 
From these figures it is clear that the profile of the line C IV 
is very similar to the spectrum of low-ionization BAL QSO 
(and different to the spectra of non-BAL QSOs). 

In general, the UV and optical spectral features of IRAS 
21219-1757 are similar to those observed in nearby low- 
ionization BAL QSOs: like IRAS 07598+6508 and PG 1700 
+518. These objects are also extreme IR + Fe n emitters 
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(Lipari 1994; Wampler 1985; Lanzeta et al. 1993; Hines & 
Wills 1995). 

In the IUE UV spectrum of IRAS 21219-1757 the follow- 
ing features were also observed: 

• The absorption feature in the C IVA1549 line ex- 
tends from ~1719 to ~1665 A, implying outflow/ejection 
velocities of ~1500 to ~10400 km s -1 . 

The emission line of C ivA1549 is also weak, as has 
been observed in BAL-QSO (Weymann et al. 1991). 

• Figs. 12a,b also show the presence of weak low- and 
high-ionization "emission" lines. In particular, Lya and 
N VA1240 are weak. 

• There are also weak emissions from O IA1303 and 
Si IV+O IV A1400. 

We note that the archive HST STIS spectra of IRAS 
21219-1757 does not cover the region where we detected the 
BAL in this IR QSO (the line C ivA1549). 

(ii) In order to compare this BAL IR QSO, with the previ- 
ous observed BAL IR QSOs, we have also included in Fig. 13 
the UV spectra of IRAS 07598+6508, PG 1700+518, IRAS 
14026+4341. 

Furthermore, it is important to note that the UV-BAL 
features observed in these IR + GW/OF + Fc II QSOs 
are very similar to those observed at high redshift in low- 
ionization BAL-QSOs (at z > 6.0; see Maiolino et al. 2003). 

(iii) In Fig. 13d the IUE spectrum of IRAS 18508-7815 is 
presented. In this extreme IR + GW/OF + Fe n QSO there 
is only evidence of a possible narrow absorption system in 
the line C IV (and close to the systemic velocity). 

(iv) For PHL 1092 (another extreme IR + GW/OF + 
Fe n QSO) we have studied mainly the region of the Mg 
n+UV-Fe II. Fig. 13e does not show evidence of Mg II BAL 
systems, in PHL 1092. However, the fit of the emission line 
in this area required to include an UV Fe II template (we 
used the UV Fe II emission of I Zw 1, from the HST spectra) 
and also a blue OF component. For this blue component a 
OF value of -8565 km s _1 was measured. 

At the optical/red wavelengths, this OF (with very high 
velocity) was not detected at the line Ha. However, we found 
3 blue OF componets with relativelly low velocities. For 
these blue components the following OF values were mea- 
sured, OF1 = -1500, OF2 = -2500 and OF3 = -3100 km 
s _1 (see Fig. 13f). These velocities are almost the same to 
those obtained for the OFs detected in H/3 (using the offset 
method; Lipari et al. 2004d: their Table 8). 

It is important to note that Dietrich et al. (2002) and 
Barth et al. (2003) already discussed that in order to obtain 
a good fit of the UV lines Mg II + Fe II in very high redshift 
QSOs, they need to include a strong blushifted component 
(they explain that this component was used without a phys- 
ical explanation). In particular, Barth et al. (2003) in their 
Fig. 2 show this strong blue OF component in the line Mg 
ii, for the Fe n-QSO SDSS J114816.64+525150.3. This ob- 
ject is one of the younger QSO known, with a redshift z = 
6.4. They measured for the Mg II OF component -of this 
distant/young QSO- a value of ~ -5700 km s _1 (Barth et 
al. 2003: their Fig. 2), which is almost 1/3 smaller than the 
value obtained for PHL 1092. 

Therefore, we strongly suggest that this type of blue com- 
ponent observed in the Mg II emission line -in very high red- 
shift QSOs- is associated with extreme OF processes. Again, 



the results obtained from QSOs at very high redshift (by Di- 
etrich et al. 2002 and Barth et al. 2003) are very similar to 
our results for IR + GW + Fe II QSOs at low redshift (where 
we detected blue OF components in Ha, H/3, [O III], [O n] 
and Mg II emission lines). However, these OF emissions are 
probably associated with very different OF processes, simi- 
lar to those found in this paper for Mrk 231. 
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4 DISCUSSION 

In this section we discuss the properties of Mrk 231 as: (i) an 
evolving elliptical galaxy; (ii) a merger with composite nu- 
clear energy; (iii) an EVOF QSO, with also a composite na- 
ture; (iv) as an exploding QSO. In addition, we analyse the 
IR colours diagram for a large sample of IR mergers/QSOs 
with galactic winds (including BAL QSOs and Mrk 231). 

4.1 The merger Mrk 231 as an evolving elliptical 
galaxy with galactic wind 

Hamilton & Kell (1987) already found that the optical R 
broad band surface brightness profiles of Mrk 231 follow the 
r 1 / 4 law (in the range 1.0 kpc < r < 10 kpc). Lipari et 
al. (1994) from high resolution NOT V broad band images 
found that the nucleus is compact at scales of 0'.'7 (~500 
pc). Soifcr et al. (2000) using Keck mid-IR 7.9 to 19.7 ^m 
broad band images found that the nucleus is compact and 
unresolved at 0'.'13 (100 pc). At radio wavelengths, Con- 
don et al. (1991) also detected an unresolved radio core, 
smaller than 0'.'25 (at 8.44 GHz), in Mrk 231. HST ACS/U, 
WFPC2/I and NICMOS/H broad band images (see section 
3 and Quillen et al. 2001) confirmed that at scales better 
than O'.'l (~80 pc, FWHM), the bright compact nuclear 
peak is unresolved. These results suggest that the nuclei of 
the original colliding galaxies have coalesced into a common 
nucleus with an AGN, and that the merger is in a very ad- 
vanced phase (a relaxed system), probably evolving to an 
elliptical galaxy. 

Therefore Mrk 231 is another candidates for a proto- 
elliptical galaxy. It is important to note that from the the- 
oretical point of view there is a wealth of literature propos- 
ing a relevant role for galactic wind and OF in the forma- 
tion and evolution of elliptical galaxies. We can highlight 
that: (i) Mathews & Baker (1971) suggested that a galac- 
tic wind can explain the deficiency of observed gas as com- 
pared to that returned by stars in ellipticals; (ii) Larson 
(1974), Vader (1986) and Kauffmann & Chariot (1998) pro- 
posed for ellipticals a galactic wind associated with the early 
star formation episode to explain their colour-magnitude 
and mass-met allicity relation; (iii) Larson (1974), Vader 
(1987) and Dekel & Silk (1986) proposed that dwarf ellip- 
ticals evolved from initially more massive elliptical galaxies 
that suffered substantial gas loss through supernova-driven 
galactic winds. 

The results obtained for Mrk 231 (and previous studies 
of mergers with OF) suggest that extreme starbursts plus 
AGN with galactic winds play an important role in galaxy 
evolution, e.g. giant expanding shells in composite galac- 
tic winds could produce BAL systems, massive starbursts 
could generate AGNs/QSOs, successive extreme starburst 
processes could transform the observed high gas-density (in 
the ISM of mergers) to high stellar-density, etc. Further- 
more, these extreme starburst+AGN with galactic winds 
-induced by mergers- could be a nearby analogy for pro- 
cesses that occurred at high redshift, when the galaxies and 
QSOs formed, i.e. similar to the first massive star forma- 
tion episodes (population III stars) in hierarchical merg- 
ers and/or primordial collapses (Larson 2003, 1999, 1998; 
Bromm & Loeb 2003; Bromm, Coppi, & Larson 1999; Scan- 
napieco & Broadhurst 2001). 



4.2 The composite nucleus of Mrk 231: an 

AGN/QSO embedded in an extreme starburst 

In order to study the nature and origin of the complex nu- 
clear outflow process detected in Mrk 231, it is important 
to analyse -previously- the composite nature of the nuclear 
energy. 



4.2.1 The AGN/QSO 

There is clear and interesting evidence of the presence of an 
AGN/QSO in the nucleus of Mrk 231. Specifically, at radio 
wavelength Preuss & Fosbury (1983) already detected a very 
compact nuclear radio source, at scale <1 pc (using VLBI 
interferometric observations) . Recently, using the Very Long 
Baseline Array (VLBA) Ulvestad et al. (1999a,b) found a 
parsec scale radio jet, with an extension of ~2.5 pc. 

At X-ray wavelength, using ASCSA data Turner (1999) 
found that the optical-to-X ray spectrum of Mrk 231 is typ- 
ical of a QSO. In addition, using Chandra observations Gal- 
lagher et al. (2002) detected that the bulk of the X-ray lu- 
minosity is emitted from an unresolved nuclear point source 
and the spectrum is remarkably hard. 

In section 3, we found clear evidence that -at least- the 
Na ID BAL system I is associated with the presence of the 
AGN+jet. 



4-2.2 The extreme starburst 

Recently, the star formation in the nucleus of Mrk 231 was 
studied at very high resolution by Bryant & Scoville (1996) 
and Downes & Solomon (1998), observing the CO gas. In 
particular, Bryant & Scoville (1996) found that the CO dis- 
tribution is elongated with kinematics properties consistent 
with a disk oriented in the east-west direction, and contain- 
ing a mass of 3 x 10 9 Mq of molecular gas (within a diameter 
of l'.'O, 814 pc). Downes & Solomon (1998) using CO inter- 
ferometer data also detected a l'/2 diameter inner disk plus 
a 3'.'0 diameter outer disk. They found extreme star forma- 
tion processes in the nuclear CO disks or rings of Mrk 231. 
Confirming previous results that suggested the presence of a 
strong nuclear starburst in this IR merger (Hamilton & Keel 
1987; Hutching & Neff 1987, and others). In particular Li- 
pari et al. (1994) also suggested the presence of an extreme 
nuclear starburst, as the main source of the strong galactic 
wind, detected in Mrk 231. These works suggested that at 
least 1/3 of the total source of the nuclear energy in Mrk 
231 is generated by the extreme starburst. 

On the other hand, continuum observations at 21 cm 
with VLBA (Carrilli, Wrobel & Ulvestad 1998) show emis- 
sion from an extended disk of 370 pc, that exhibits absorp- 
tion against the radio emission from the inner disk (r < 
160 pc). Further radio continuum observations reported by 
Taylor et al. (1999) trace this emission out to l'.'O. They in- 
terpret this disk radio emission also as due to strong star for- 
mation activity with a SFR of 60-200 M yr _1 . At near-IR 
wavelengths, the nuclear spectrum of this galaxy shows an 
anomalously large Paa/Br7 ratio (Cutri et al. 1984), which 
imply that the electron density in the broad line clouds of 
Mrk 231 is higher than the standard values of the BLR in 
AGNs/QSOs. 
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In section 3.4, it was found that extended shocks (as- 
sociated with multiple expanding bubbles) are also an im- 
portant source of energy and ionization, in the nuclear and 
circumnuclear regions. In addition, we found evidence of a 
multiple explosive event in Mrk 231, which is in part associ- 
ated with the extreme starburst, probably generated during 
the merger process. 

At least two strong starbursts were generated in the 
merger process of Mrk 231: one associated with the origin of 
an extended post-starburst population (associated with the 
underlain hot stellar population in the optical spectrum), 
plus a young and dusty nuclear extreme starburst. 



4.3 The composite nuclear outflow + BAL 
systems in Mrk 231: an exploding QSO? 

In section 3 very interesting results, about the OF and BAL 
systems in Mrk 231 were found. Specifically, it was detected: 
(i) multiple concentric expanding bubbles, (ii) a blue Ha- 
cmission bump at the ejection velocity of the BAL I, and 
(iii) a light curve for the BAL-III system showing a steep 
increase and exponential fall (very similar to the shape of a 
SN light curve). 

In order to analyse these new results, it is important to 
remark previous studies at radio wavelengths, which are rel- 
evant for the discussion of the OF in Mrk 231. In particular: 

(i) VLBI images at 1.7 Ghz (Neff & Ulvestad 1988) and 
VLBA images at 2.3 Ghz (Ulvestad et al. 1999a) show a 
north-south triple structure with a central unresolved core 
and two symmetric resolved lobes; with a total extension of 
^40 pc. This radio structure is elongated (at PA ~00°) in 
the direction perpendicular to a 350 pc HI + CO starburst 
disk/ring. 

(ii) VLA images at 4.9 Ghz (Baum et al. 1993) and at 1.5 
Ghz (Ulvestad et al. 1999a) show a very large structure of 
~35-50 kpc, also elongated in the direction north-south. 

(iii) Using the VLBA at 15.3 GHz, Ulvestad et al. 
(1999a, b) found a parsec scale radio jet, with an extension 
of ~2.5 pc, low apparent speed, and at the position angle 
PA ~65°. 



4-3.1 Multiple expanding shells 

The presence of multiple concentric expanding bub- 
bles/shells, with centre in the nucleus and with highly sym- 
metric circular shape could be associated mainly with gi- 
ant symmetric explosive events. These giant explosive events 
could be explained in a composite scenario: where mainly the 
interaction between the starburst and the AGN could gen- 
erate giant explosive events. In particular, Artymowicz, Lin, 
& Wampler (1993) and Collin & Zahan (1999) already anal- 
ysed the evolution of the star formation (SF) close to super 
massive black hole (SMBH) and inside of accretion disks. 
They suggested that the condition of the SF close to the 
AGNs could be similar to those of the early/first SF events, 
where it is expected giant explosive process generated by 
hyper-novae (with very massive progenitors: M ~100-200 
M ). 

For the shell SI, Lipari et al. (1994) already proposed 
that the morphology and kinematics of the blue arc (located 



to the south of the nucleus) are in agreement with the re- 
sults of glactic winds hydrodynamic models. Specifically, the 
starburst GW models of Suchkov et al. (1994) predict after 
8 Myr in the blowout phase: blue arcs with radius r ~ 3 kpc. 
We noted that after 8 Myr the starburst is in the type II SN 
phase. 

It is important to remark that exactly in the direction 
of this arc (north-south from the nucleus, PA ~00°) Neff &: 
Ulvestad (1988) and Ulvestad et al. (1999a) found a elon- 
gate structure of ~40 pc, which is perpendicular to the CO 
starburst disk/ring. Furthermore, Baum et al. (1993) and 
Ulvestad et al. (1999a) detected a large structure (of ~35- 
50 kpc), also elongated in the direction north-south. These 
results are all clearly consistent with the blowout or rupture 
phase of the bubble SI. 

Following our study of the superbubble in NGC 5514 
(Lipari et al. 2004d) the dynamical timescale of the more 
extended shells SI and S2 in Mrk 231 were derived, using 
the relation: 

tdyn = 10 XlO Rbubble,kpc V bubble , 1000 ^ r 

where Rbubbie,k P c and V^ 1 We . il000 are the linear dimen- 
sion of the bubble and the velocity of the entrained material 
(in units of kpc and 1000 km s _1 , respectively). The follow- 
ing values were obtained for the shells SI and S2: 

• Uyn-si ~ 4.8 x 10 6 yr; 

• td yn -S2 ~ 3.0 x 10 6 yr. 

Then a period of ~ 1.8 x 10 6 yr is the separation be- 
tween the two explosive events that probably generate the 
shells SI and S2. 

From the 2D spectra a mean value of Ha flux of 0.5 
x 10~ 13 erg cm~ 2 s _1 was obtained, for the OF associated 
with the shell SI. Using the relations given by Mendes de 
Olivera et al. (1998) and Colina et al. (1991) a value for 
the ionized-gas mass, in the OF associated with the shell 
SI, of 2.0 x 10 6 M Q was found. Using this value of ionized- 
gas mass, we derived for the kinetic energy of the OF-S1, 
Ekinof-si ~ 0.5 x Mof-ig x (Vof-si) 2 = 2.0 x 10 54 erg. 
This high level of energy, obviously required the presence of 
multiple SNe events, or an unusual type of "giant SN" or 
hypernovae. Heiles (1979) already suggested that this last 
alternative need to be considered seriously, in order to ex- 
plain supershells with enegy > 3.0 x 10 52 erg (even if these 
energies of supershells are hundred of times larger than that 
available from a single and "standard" SNe) . 

In section 4.3.4 the origin of this multiple shells will be 
discussed in detail. 



4-3.2 The blue Ha-emission bump at velocity of BAL I 
system 

In section 3 another interesting results was found: a blue 
Ha-emission bump at the ejection velocity of the BAL I 
(and in a area located very close to the nucleus and at the 
same PA of the radio jet). Furthermore, in the same di- 
rection where we detected this bump (at PA ~60-70°) at 
l'.'3-2'.'2, we also found multiple narrow emission line compo- 
nents, with "greatly" enhanced [Nil] /Ha ratio (very similar 
to the spectra of jets bow shocks). These results are consis- 
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tent with a scenario where the BAL-I system is generated 
in OF clouds associated with the parsec scale jet. 

It is important to note, that de Kool et al. (2001, 2002) 
found a very interesting result -and physically similar to 
those found for Mrk 231- in their Keck high resolution 
spectroscopic study of AAL and BAL, in the QSOs FIRST 
J104459.6+365605 and FBQS 0840+3633. The absorption 
lines of these QSOs cover a ranges of velocities of -200 to - 
1200, -3400 to -5200 km s" 1 ; and the width of the individual 
absorption lines ranges from 50 to 1000 km s _1 (these arc 
values similar to those measured in Mrk 231). They found 
that the distances between the AGN and the region where 
the OF gas generate the AAL and BAL line are ~700 and 
~230 pc, respectively. Therefore, for Mrk 231 (BAL system- 
I), J104459.6+365605 and FBQS 0840+3633 the distances 
found between BAL or AAL forming regions and the contin- 
uum source (AGN) are large, of: ~ 200-1000 pc (BALs and 
AALs are generally thought to be formed in OF at a much 
smaller distance from the nucleus; de Kool et al. 2001). 

4-3.3 The origin of the variability of the BAL III system 

Boroson et al. (1991) and Kollatschny et al. (1992) already 
explained some physical mechanism that could explain vari- 
ations in the absorption line strengths. Here we re-analyse 
these and new mechanisms, according to the results obtained 
in section 3. 

The principal mechanisms that can explain the vari- 
ability of the NA ID BAL III system, can be enumerate as 
follow: 

(i) Changes in the intensity of the continuum source (and 
the ionizing parameter in the absorbing region). Since the 
continuum flux remained constant in the period in which 
the system III varied, the first possible explanation can be 
excluded. 

(ii) Changes in the ionizing parameters and density in 
the absorbing clouds when they move away from the centre. 
This second possibility have been excluded by Boroson et 
al. (1991), since the relative change in distance can not ex- 
plain the strong change and the systematic increasing and 
decreasing epochs. 

(iii) Motion of absorption clouds with traverse velocity 
occulting the continuum source. This was in the past the 
best explanation for part of the light curve of BAL III, 
observed between 1980 and 1992. However, the new obser- 
vation -presented in section 3.5- exclude this explanation, 
sincethe LC fall is clearly asymmetric/exponential. 

(iv) An explosive scenario for the origin of the BAL III 
system could explain: the shape of the light curve variabil- 
ity, and also the presence of multiple concentric expanding 
bubbles/shells (with circular shape). In the next section we 
will discuss in detail this explosive scenario and the probable 
origin. 

4-3.4 The exploding scenario for Mrk 231 and BAL 
systems 

In section 3 strong evidence for a composite origin of the OF 
in the nuclear region of Mrk 231 were detected. Furthermore, 
for IR mergers/QSOs with EVOF a similar result was found 
(see for details section 1 and Lipari et al. 2004a,d). Thus 



it is probably that also very energetic explosive events are 
originated in the interaction between extreme starbursts and 
AGNs. We already explained that Artymowicz et al. (1993) 
and Collin & Zahan (1999) suggested that the conditions of 
the SF close to the AGNs+accretion disks could be similar 
to those of the early/first SF events, where it is expected 
giant explosive processes generated by giant SNe. 

On the other hand, different works also proposed that 
the BLR could be associated with OF processes in accre- 
tion disks, extended stellar envelopes, etc (see for refer- 
ences/review Sulentic, Marziani, & Dultzin-Hacyan 2000; 
Scoville & Norman 1996; Dyson & Perry 1992; Terlevich et 
al. 1992). Thus an interesting questions is: could an explo- 
sive event generate a very unusual spectrum similar to that 
detected in the nuclear region of Mrk 231?. Figure 14 shows 
the superposition of the spectrum of the unusual radio SN 
of type II 1979c (observed in 1979 June 26.18; Branch et al. 
1981). Only using colours we can distinguish each spectrum, 
since they are almost identical. Thus, an more constant OF 
(than a single/standard SN) could explain even the spec- 
trum of the BLR in Mrk 231. 

In the composite (starburst+AGN) scenarios, two main 
theoretical models for the origin of BAL systems were pro- 
posed: (i) for IR dusty QSOs/galaxies, in the outflowing 
gas+dust material the presence of discrete trails of debris 
(shed by individual mass-loss stars) produce the BAL fea- 
tures (Scoville & Norman 1996); and (ii) in SN ejecta, which 
are shock heated when a fast forward shock moves out into 
the ISM (with a velocity roughly equal to the ejecta) and a 
reverse shock accelerates back and moves towards the explo- 
sion centre; the suppression of red-shifted absorption lines 
arise since SN debris moving toward the central source are 
slowed down much more rapidly -by the wind- than is mate- 
rial moving away (Perry & Dyson 1992; Perry 1992). The 
presence of large galactic-scale OF, shells or rings in IR 
QSOs could be a third explanation -in a composite sccnario- 
for the origin of BAL systems, in these objects. 

It is interesting to recall that in their study of PG 
1700+518, Hazard et al. (1984) already suggested that: high 
redshift low-ionization BAL-QSOs could be explained by a 
violent ejection during the first onset of the QSO activity, 
similar to a "giant SN explosion". This approach is very sim- 
ilar to that proposed in this paper for Mrk 231 (where "giant 
expandig superbubbles" were foundd). Therefore, it is impor- 
tant to study in detail the contribution of extended outflow 
associated with superwind/superbubble regions in order to 
explain the origin of some BAL systems and even the BLR of 
Mrk 231 (see Lipari 1994; Lipari et al. 1994; Terlevich et al. 
1992; Dyson, Perry & Williams 1992; Voit, Weymann & Ko- 
rista 1993). In general, the results and conclusions presented 
in Section 3 appear to be in good agreement with this un- 
usual explosive scenario proposed in this paper for Mrk 231: 
where the strong IR+Fe II emission and BAL system could 
be linked -in part- to a violent extreme-starburst+AGN. 

Finally, the serendipitous finding that a high number of 
QSOs show BAL systems at low redshift (in IR + Fe n emit- 
ters) and at very high redshift (in Fe II emitters; Maiolino 
et al. 2003, 2004) is very interesting. On the other hand, the 
relation between BAL, IR, GW and Fe II emission is not 
well understood. Fe II emission (with low ionization level) 
could be originate in warm regions, and thus obscured from 
the direct ionizing UV photons. In the case of strong Fe II 
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emitters, the obscuring material may be in the form of an ex- 
panding shell. If so, BAL+Fe Il+IR objects are not ordinary 
AGN, rather, they may be associated with early starburst 
activities. If some absorptions are the result of outbursts, 
the mass involved should be much larger than a typical su- 
pernova, since the Fe II emission and "some" BAL systems 
do not fall in years. However, in order to understand this 
type of giant outbursts (i.e., from hypernova, SN explosion 
from population III stars, SN in accretion discks of AGNs, 
etc) it is required more detailed theoretical studies. The re- 
sults found in this paper for Mrk 231 (and similar IR QSOs) 
support the reality of these giant explosive processes. 



4.4 Mrk 231 and extreme IR + GW/OF + Fe II 
+ BAL QSOs as transition/young QSOs 

In this section, we will analyse the properties of Mrk 231 
together with the other members of the interesting group of 
IR + GW/OF + Fen + BAL transition/young QSOs. 



4-4.1 The IR colours diagram for a small sample of IR + 
GW/OF + Fell + BAL QSOs 

In general, the IRAS colour-colour diagrams have been used 
as an important tool to detect and discriminate different 
types of activity in the nuclear/circumnuclear regions of 
galaxies (e.g., Seyfert and starburst activity; see de Gripj 
ct al. 1985, 1987; Sekiguchi 1987; Rowan-Robinson & Craw- 
ford 1989). Lipari (1994) already found that the IR colours 
(i.e., IR energy distribution) of ~10 extreme IR + Fe 11 QSOs 
(e.g., Mrk 231, IRAS 07598+6508, IRAS 17002+5153, IRAS 
14026+4341, IRAS 00275-2859, I Zw 1, IRAS 18508-7815, 
IRAS 21219-1757, Mrk 507, Mrk 957, etc) are distributed 
between the power law (PL) and the black-body (BB) re- 
gions: i.e., the transition area. On the other hand, the low 
and strong IR-Fe II emitters are located mainly in the PL 
region. 

In particular, we detected that Mrk 231 and 
IRAS 07598+6508 (the nearest IR + GW/OF + Fe 11 + BAL 
QSOs) have a close position in this diagram: near to the BB 
area; thus showing both systems strong starburst compo- 
nents. Recently, Canalizo & Stocton (2001) confirmed that 
the host galaxies of both QSOs have strong starburst popu- 
lations (using Keck spectroscopy). In general, the extreme Fe 
II + IR emitters show weaker "Seyfert/ AGN" components 
than the low Fe II QSOs. 

Furthermore, we found that the IR properties/colours 
of narrow-line Seyfert 1 IR AGN with extreme Fe II emission 
(like Mrk 507 and Mrk 957), lie also in the transition area 
between the BB and PL regions, and close to standard disk 
galaxies (Rowan-Robinson & Crawford 1989) and starburst 
objects (Sekiguchi 1987). 

It is important to remark that of a total of ~10 IR 
transition objects of this original sample, the first 4 systems 
are BAL IR QSOs. Therefore, we already suggested that 
BALs IR QSOs (like Mrk 231, IRAS 07598+6508, IRAS 
17002+5153 and IRAS 14026+4341) could be associated 
with the young phase of the QSO activity. 



4-4-8 The IR colours diagram for a large sample of IR 
mergers and IR QSOs with galactic winds 

In this paper, using our data base of more than 50 IR Merg- 
ers and QSOs with galactic winds and using for compar- 
ison the large sample of standard PG QSO (from Boro- 
son & Green 1992) we will expand our previous study. 
Figure 15 and Table 4 show the IR energy distribution 
[spectral indexes q(60, 25) vs. a(100, 60); where a(A2, Al) 
= -log[F(A2)/F(Al)]/log[A2/Al]] for: (i) IR mergers and IR 
QSO with GW (originally 51 IR systems, from Table 4); (ii) 
standard QSOs from the PG QSOs sample of Boroson & 
Green (1992; originally 87 PG QSOs that have z <0.5). 

IR fluxes-densities, in the bands of 12, 25, 60 and 100 
/jm, were obtained from the IRAS and ISO Archival Cat- 
alogue (using NED). Only objects with a good detection 
in the three required bands have been included. Also, the 
localisation of the three main regions in this colour-colour 
diagram (i.e., the QSOs/Seyferts, Starbursts, and powerful 
IR galaxies areas) have been plotted. An inspection of this 
diagram clearly shows the following: 

(i) All the IR mergers with LVOF are located very close 
to the BB and starburst area. 

(ii) Almost all the IR QSOs with EVOF are located in 
the transition region. 

The only object, NGC 3079, not following this trend is 
the only spiral galaxy in this sample. 

(iii) The standard QSOs and radio QSOs are located 
around the PL region. 

(iv) All the BAL IR QSOs are located in the transition 
region, in almost a clear sequence: from Mrk 231 (close to 
the BB area) -> IRAS 07598+6508 -> IRAS 21219-1757 
IRAS/PG 17072+5153 and IRAS 14026+4341 (close to the 
PL area) — > standard QSOs. 

Therefore, these results confirm our previous finding 
(obtained from a small sample of IR galaxies): in the sense 
that IR QSOs are probably "transition" objects, between IR 
mergers and standard QSOs. However, it is important to re- 
mark that it is required more detailed studies of individual 
IR systems (like Mrk 231) in order to determine the dif- 
ferent paths of possible evolution (see Farrah et al. 2001; 
Rigopoulow et al. 1999). 



44.3 The possible relation between IR + GW/OF + Fe 11 
+ BAL transition/young QSOs and very high 
redshift BAL QSOs 

It has been proposed that extreme starburst + galactic 
wind processes associated with IR mergers could play a rel- 
evant role in the formation and evolution of galaxies and 
QSOs/AGNs, i.e. in their structure, kinematics, metallicity, 
etc. (see sections 1 for references). Furthermore, recent de- 
tailed observations and theoretical studies have confirmed 
that OF, galactic winds, BAL, large amount of gas+dust 
and strong Fe 11 emission are important components and 
processes at high redshift (z ~ 4-6), when the galaxies and 
QSOs formed (Frye, Broadhurst, & Benitez 2002; Ajiki et 
al. 2002; Pettini et al. 2001; Dawson et al. 2002; Carilli et 
al. 2004a,b; Solomon et al. 2003; Taniguchi & Sioya 2000; 
Barth et al. 2003; Iwamuro et al. 2002; Freudling, Corbin, 
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& Korista 2003; Bertoldi et al. 2003a,b; Omont et al. 2001; 
Maiolino et al. 2004a,b). 

Specifically, cm and mm observations of thermal molec- 
ular line emission from high redshift QSOs reveal massive 
gas reservoirs (10 10 to 10 11 Mq) required to fuel extreme 
star formation (Carilli et al. 2004a,b; Solomon et al. 2003; 
Carilli & Blain 2002; Carilli, Menten & Yung 1999). In ad- 
dition, cm and mm observations also reveal in high redshift 
QSOs very large amount of dust of 10 8 M . The presence 
of this large mass of dust is associated mainly with massive 
starbursts concurrent with AGNs, which implied star for- 
mation rates greater than 10 3 M yr _1 (Carilli et al. 2000; 
Bertoldi et al. 2003a,b; Omont et al. 2001; Cox et al. 2002). 

On the other hand, Maiolino et al. (2004a,b, 2003) pre- 
sented near-IR spectra of eight of the more distant QSO 
(at 4.9 < z < 6.4). Half of these QSOs are characterised 
by strong UV BAL systems (at C IV, Mg II, Si IV, Al III 
lines). Although the sample is small, the large fraction of 
BAL QSOs suggest that the accretion of gas, the amount 
of dust and the presence of OF process are more larger (in 
these objects) than in standard QSO at z < 4.0. They also 
suggested that the very high amount of dust was generated 
by early explosions of SNe (Maiolino et al. 2004b). 

Finally, it is important to remark the similar properties 
found in IR + GW/OF + Fc II + BAL QSOs at low redshift 
(like Mrk 231) and very high redshift BAL QSOs (at z ~ 
6.0; Maiolino et al. 2003, 2004a,b; Carilli et al. 2004a,b). 
According to these similarities, we propose that the phase 
of young QSO could be associated with the following main 
processes: 

(i) In young QSOs with extremely large amount of gas 
(concentrate in their nuclear region), the accretion rate of 
gas -by the SMBH- could be extremely high (see Maiolino 
et al. 2004a). 

(ii) In addition this extremely large amount of molecular 
gas could generate extreme starbursts; and the presence of 
AGNs could increase the SF close to the nucleus. Specially 
in the accretion disks, with properties -of the SF- similar 
to the population III of stars. In these extreme starbursts 
-close associated with QSOs/ AGNs- it is expected giant SN 
explosions. 

(iii) In young and distant QSOs the very high number 
of BAL detections suggest that composite OFs (or EVOFs) 
play a main role in their evolution. 

Furthermore, this scenario is in close agreement with 
the results of the study of proto-QSO evolution and black 
hole growth (using hydrodynamic models). In particular, 
Kawakatu et al. (2003) found that a ULIRG phase (where 
the host galaxy is the dominant source of luminosity: i.e. 
IR galaxies/mergers with starbursts) precedes a galactic 
wind epoch: i.e., composite young IR+OF/GW+BAL merg- 
ers/ QSOs. Which would be a transition state to the AGN- 
dominated phase: i.e., the standard QSO phase. This evolu- 
tive path is almost identical to the sequence found for BAL 
IR+Fe ii QSOs, in Figure 15. 
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5 SUMMARY AND CONCLUSIONS 

We have presented in this paper optical 2D and ID 
spectroscopy (obtained at La Palma/WHT, HST, IUE, 
ESO/NTT, KPNO, APO and CASLEO observatories) and 
deep HST broad band images of Mrk 231 and similar IR + 
GW/OF + Fe II QSOs. The main results and conclusions 
can be summarised as follows: 

(i) A detailed study of the kinematics, morphology, physi- 
cal conditions and ionization structure of the extreme galac- 
tic superwind + bubble (detected previously in the nu- 
clear and central regions of the IR merger/QSO Mrk 231) 
is presented. This study is based mainly on INTEGRAL 
two-dimensional INTEGRAL spectroscopy (obtained at La 
Palma 4.2 m William Herschel Telescope) plus HST near- 
UV/ACS, optical/WFPC2 and near-IR/NICMOS broad 
band images. 

(ii) Inside and very close to the nucleus the 2D spectra 
show the presence of an OF emission bump in the blend 
Ha+[N n], with a peak at the same velocity of the main 
BAL-I system (VEjectionBAL-i ~ -4700 kms -1 ). This bump 
was more clearly detected in the area located at 0'.'6-l'.'5 
(490-1220 pc), to the south-west of the nucleus core, showing 
a strong and broad peak. In addition, in the same direction 
(at PA ~60, i.e. close to the PA of the small scale radio 
jet) at l'.'7-2'.'5 (1340-2030 pc), we also detected multiple 
narrow emission line components, with "greatly" enhanced 
[Nil] /Ha ratio (very similar to the spectra of bow shocks 
in jets, already found in the Circinus galaxy). Thus, these 
results suggest that probably the BAL-I system is generated 
in OF clouds associated with the parsec scale jet. 

(iii) The HST images show 4 (or 5) nuclear bubbles or 
shells with radius r ~ 2.9, 1.5, 1.0, 0.6 and 0.2 kpc (3'.'5, 
l'/8, l'/2, 0'.'7 and 073). For these bubbles, the 2D Ha veloc- 
ity field (VF) map and 2D spectra show: (i) At the border 
of the more extended bubble or shell (SI), a clear expansion 
of the shell with blueshifted velocities (with circular shape 
and at a radius r ~ 3'.'0). This bubble shows a rupture arc 
-to the south- suggesting that the bubble is in the blowout 
phase. The axis of this rupture or ejection (at PA ~00°) is 
coincident with the axis of the intermediate and large scale 
OF detected at radio wavelengths (r ~ 0'.'6 ~ 50.0 pc at 2.3 
GHz, and r ~ 307 ~ 25.0 kpc at 1.5 GHz), (ii) In addition 
in the 3 more external bubbles (SI, S2, S3), the 2D WHT 
spectra show multiple emission line components with OF 
velocities, of (VoFBubbic-si) = [—(410 — 650) ± 30]kms _1 , 
(V FBubbie-S2) = [-500 ± 30] km s _1 , and (VoPBubbie-S3> = 
[—230 ±30] kms" 1 . (iii) In all the circumnuclear region (l'/8 
< r < 5"), the [Nil] /Ha and [S ll]/Ha narrow emission line 
ratios show high values (> 0.8), which are consistent with 
LINER/OF processes associated with fast velocity shocks. 
Therefore, we suggest that these giant bubbles are associ- 
ated with the large scale nuclear OF component, which is 
generated -at least in part- by the extreme nuclear star- 
burst: i.e., type II SN explosions. 

(iv) High resolution HST ACS, WFPC2 and NICMOS 
broad band images (using the filters F330W, F439W ~B, 
F814W ~J, and F160W ~H) were combined with the 2D 
WHT spectra to study the physical properties of the mul- 
tiple bubbles/shells. Specially, we study the properties of 
several blue/young star forming regions, located mainly at 
the border of these shells. 



(v) Finally, the variability of the short lived BAL-III 
Na ID system was studied, covering almost all the period 
in which this system appeared (between ~ 1984-2004). We 
found that the BAL-III light curve (LC) is clearly asymmet- 
ric with: a steep increase, a clear maximum and an exponen- 
tial fall (similar to the shape of a SN LC). Previously the 
nature of this LC was discussed assuming a symmetric or 
Gaussian form. In this paper we discuss the origin of this 
BAL-III system, mainly in the frame work of an extreme 
explosive event (at parsec scale), probably associated with 
giant SNe or hyper-nova, with very massive star progenitors 
and located close to the AGN and/or in an accretion disk. 

Using our data base for more than 50 IR Mergers and 
QSOs with galactic winds and using for comparison the large 
sample of standard PG QSO (from Boroson & Green 1992), 
we confirm our previous finding (obtained from a small sam- 
ple of IR galaxies): in the sense that IR QSOs are proba- 
bly "transition" objects, between IR mergers and standard 
QSOs. 

Finally, new observations of UV-BAL systems of IR + 
GW/OF + Fe II QSOs were analysed. This study shows a 
new BAL IR QSO and suggest /confirm that these objects 
could be nearby young BAL QSOs, similar to those detected 
recently at z ~ 6.0. We propose that the phase of young 
QSO is associated with: accretion of large amount of gas 
(by the SMBH) + extreme starbursts + extreme composite 
OFs/BALs. 
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Table 1. Journal of observations of Mrk231 and similar extreme 
IR + GW/OF + Fc II QSOs 



Object 


Date Telescope/ 


Spectral region 


Expos, time Comments 




instrument 




[s] 


Mrk 213 


1991 Feb 15 2.15 m KPNO/GoldCam 


AA3350-5200 A 


900x2 PA = 90, slit width = 1"5 




2.15 m KPNO/GoldCam 


AA5100-7100 A 


900x2 



Mrk 231 


1991 May 11 


2.5 m NOT 


V 


1200x3 


(FWHM) 


= 0"7 


Mrk 231 
Mrk 231 
Mrk 231 


1992 Nov 27 
1992 Nov 27 
1996 Nov 21 


HST/FOS 
HST/FOS 
HST/FOS 


G190H, AA1275-2320 A 
G270H, AA2225-3295 A 
G160L, AA1150-2300 A 


5760 
2880 
770 


archival 




Mrk 231 
Mrk 231 


1995 Oct 23 
1995 Oct 23 


HST/WFPC2 
HST /WFPC2 


F439W, AA4283/464 A 
F814W, AA8203/1758 A 


2226 
712 


(FWHM) 


= O'.'l, archive 


Mrk 231 


2003 Mar 17 


HST/ACS 


F330W, AA3354/588 A 


1140 


(FWHM) 


= O'.'l, archival, NIC2 



Mrlr 9^1 
IVIIK ZOl 


iijuo ocp uy 


no l 1 IN loivl KJiJ 


r 1UU VV , AA1.UU/U.4:U flLLL 


riZLfl 


(FWHM) = 0'.'22, archival, 


121219-1757 


1989 Jun 29 


2.15m CASLEO/Z-machine 


AA4700-7200 A 


1500x3 


aperture: 3"x6" 




121219-1757 


1993 Jul 13 


2.15m CASLEO/UCS 


AA6500-7000 A 


1800x2 


PA 90° slit width= 


2'.'0 


TOI Ol O 1 1K1 

izi.ziy-1 tot 


iyy ( iviar lo 


O ic™ r A CT T?H /TTPO 
Z.LOIH L^AoL/rL/U/ ULo 


\ \/iooo 7t;nri A 
AA4UUU- (OUU A 


1 COO \y O 

loUUX z 


rJ\ yU Silt W Kit 11= 


1. O 


121219-1757 


1995 Oct 24 


2.15 m KPNO/GoldCam 


AA4500-9000 A 


1800x2 


FA = yu, slit width 


= 1.0 


121219-1757 


1995 May 26 


IUE/LWP30773 


AA1150-1980 A 


9000 






121219-1757 


1995 May 26 


IUE/SWP54757 


AA1900-3290 A 


112000 






121219-1757 


1999 Jun 03 


HST/STIS 


G160L, AA1150-1550 A 


2420 


cxL L-Ill VL 




107598+6508 


1991 Feb 15 


2.15 m KPNO/GoldCam 


AA3350-5200 A 


1800x2 


PA = 90, slit width 


= l'/5 


107598+6508 


1991 Feb 15 


2.15 m KPNO/GoldCam 


AA5100-7100 A 


1800x2 






107598+6508 


1998 Dec 10 


HST/STIS 


G140L, AA1000-1400 A 


2570 


archive 




107598+6508 




IUE/SWP37199 


AA1150-1980 A 




archive 




107598+6508 




IUE/SWP40376 


AA1150-1980 A 




archive 




107598+6508 




IUE/LWP16445 


AA1900-3290 A 




archive 




107598+6508 




IUE/LWP19437 


AA1900-3290 A 




archive 




114026+4341 


1990 Sep 19 


2.15 m KPNO/GoldCam 


AA4350-5700 A 


1800x2 


PA = 90, slit width 


= 1'.'5 


114026+4341 


1994 Jul 17 


HST/FOS 


G190H, AA1570-2300 A 


978 


archive 




114026+4341 




IUE/LWP20600 


AA1150-1980 A 




archive 




114026+4341 




IUE/SWP36278 


AA1900-3290 A 




archive 




117002+5153 


1990 Sep 19 


2.15 m KPNO/GoldCam 


AA4350-5700 A 


1800x2 


PA = 90, slit width 


= l'/5 


117002+5153 


1993 Mar 10 


HST/FOS 


G190H, AA1600-2300 A 


2026 


archive 




117002+5153 




IUE/LWP4596 


AA1150-1980 A 




archive 




117002+5153 




IUE/LWP6655 


AA1150-1980 A 




archive 




PHL1092 


1993 Jul 15 


3.5 m ESO NTT 


AA3200-5000 A 


1500x3 


(FWHM) = l'.'O 




PHL1092 


2004 Sep 08 


3.5 m APO 


AA3500-5600 A 


1800x2 


(FWHM) = l'.'O 




PHL1092 


2004 Sep 08 


3.5 m APO 


AA5300-10000 A 


1800x2 






118508-7815 


1991 Jul 15 


3.6 m ESO/EFOSC 


AA3600-5590 A 


1800x2 


PA = 90, slit width 


= 1'.'5 


118508-7815 


1991 Jul 15 


3.6 m ESO/EFOSC 


AA5050-7000 A 


1800x2 






118508-7815 


1991 Jul 15 


3.6 m ESO/EFOSC 


AA6600-9800 A 


1800x2 






118508-7815 


1991 Jul 18 


IUE/LWP20841 


AA1150-1980 A 


25200 






118508-7815 


1995 May 27 


IUE/SWP54772 


AA1900-3290 A 


24300 
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Table 2. Main properties of the circumnuclear knots (associated 
with the shells) 



Knots 


a 


p ^toa vv 


TO 

r o 1 1 vv 


B-I 


Tf a 

e J J 


Comments 




L J 


L J 1 oj 


[mag] 




[pel 




1 


1.10 


1.68 22.23 


21.12 


1.11 


66 


Shell SI 


2 


-1.52 


-0.90 23.27 


20.86 


2.41 


<20 


Shell S2 


Q 
O 


1 in 

1.1U 


-U.10 Zz.oo 


on fifi 


l.OZ 




QVir.ll QQ 
OllCll OO 


4 


2.12 


-1.45 23.98 


22.60 


1.38 


42 


Shell SI 


5 


1.70 


-2.97 21.31 


20.50 


0.81 


<20 


Shell SI 


6 


1.56 


-2.74 22.97 


22.28 


0.69 


71 


Shell SI 


7 


1.01 


-3.70 22.69 


21.47 


1.22 


38 


Shell SI 


11 


-0.46 


-3.66 21.17 


20.66 


0.51 


136 


Shell SI 


12 


-0.69 


-3.61 21.46 


20.70 


0.76 


31 


Shell SI 


14 


-1.20 


-3.15 21.79 


20.33 


1.46 


40 


Shell SI 


29 


-2.30 


1.82 23.26 


21.06 


2.20 


46 


Shell SI 


Table 3. Relative variability of the Na ID BAL— III system 


Date 




Julian Date Intcns. Ratio 


a 


References 






2442507+ 










Previous Data 












04/05/1975 


0000 


0.000 


0.000 


(1) 


and (2) 


05/01/1981 


2218 


<0.005 


0.000 


(1) 


and (3) 


12/21/1984 


3548 


0.040 


0.010 


(1) 


and (4) 


02/02/1987 


4321 


0.079 


0.006 


(1) 




07/06/1988 


4841 


0.090 


0.004 


(5) 


and (1) 


01/10/1992 


6124 


0.052 


0.007 


(1) 




05/21/1992 


6257 


0.050 


0.003 


(1) 




New Data 












02/01/1991 


5781 


0.062 


0.008 


(6) 


and (7) 


02/15/1991 


5795 


0.060 


0.007 


(6) 


and (8) 


04/22/1994 


6558 


0.035 


0.003 


(6) 


and (9) 


02/28/2001 


9461 


0.011 


0.004 


(6) 


and (10) 



Notes of Table 2: 
a : From Surace et al. (1998). 

Notes of Table 3: 
References: (1) Kollatschny et al. (1992), (2) Bocksemberg 
et al. (1975), (3) Hamilton & Keel (1987), (4) Schmidt & 
Miller (1985), (5) Boroson et al. (1991), (6) This paper, (7) 
Boroson & Mayers (1992), (8) Lipari et al. (1994), (9) Forster 
et al. (1995), (10) Rupke et al. (2002). 
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Table 4. IR Mergers and IR QSOs/AGN with low and 
high/extreme velocity outflow 



Object 


Vof 1 


z 


^FIR 




Nuclear 


Morph 


RF 


BAL 


CKgO 25 


ai00,60 


OF 


OF 




km/s 




logL/L© 




activity 


Type 










Type 


Reference 


Low velocity OF 


























Arp 220 


-450 


0.0183 


12.18 


87 


L+SB 


PM 


— 


— 


-2.94 


-0.16 


EL 


Hcckman et al.1987 


IRAS00182-7112 


-450 


0.3270 


12.90 


200 


L+SB 


(M) 


— 


— 


-2.51 


+0.05 


EL 


Heckman et al.1990 


IRAS03250+1606 


-431 


0.1290 


12.06 


— 


L 


M 


— 


— 


— 


-0.49 


AL 


Rupke et al.2002 


IRAS03514+1546 


-200 


0.0210 


11.10 


— 


SB 


M 


— 


— 


-2.32 


-0.79 


AL 


Heckman et al.2000 


IRAS09039+0503 


-656 


0.1250 


12.07 


— 


L 


M 


— 


— 


— 


-0.65 


AL 


Rupke et al.2002 


IRAS11387+4116 


-511 


0.1490 


12.18 


— 


SB 


OM 


— 


— 


— 


-0.76 


AL 


55 


IRAS23128-5919 


-300 


0.0449 


12.60 


8 


L+SB 


M 


— 


— 


-2.18 


-0.03 


EL 


Li'pari et al.2000b 


Mrk 266 


-300 


0.0290 


11.37 


8 


SB 


PM 


— 


— 


-2.14 


-0.80 


EL 


Wang et al. 1997 


Mrk 273 


-600 


0.0385 


12.14 


36 


L+SB 


PM 


— 


— 


-2.58 


+0.03 


EL 


Colina et al.1999 


NGC 1614 


-400 


0.0155 


11.61 


18 


L+SB 


M 


— 


— 


-1.70 


-0.02 


EL 


Ulrich 1972 


NGC 2623 


-405 


0.0185 


11.55 


17 


L+SB 


M 


— 


— 


-2.96 


-0.37 


EL 


Li'pari et al.2004d 


NGC 3256 


-370 


0.0094 


11.57 


9 


SB 


MM 


— 


— 


-1.95 


-0.52 


EL 


Li'pari et al. 2000a 


NGC 3690 


-300 


0.0103 


11.91 


23 


SB 


PM 


— 


— 


-1.71 


-0.07 


EL 


Heckman et al.1990 


NGC 4039 


-365 


0.0056 


10.99 


6 


L+SB 


PM 


— 


— 


— 


— 


EL 


Li'pari et al.2000b 


NGC 5514 


-320 


0.0245 


10.70 


1 


L+SB 


PM 


— 


— 


-1.96 


-1.14 


EL 


Li'pari et al.2004d 


Extreme OF 


























IRAS01003-2238 


-770 


0.1180 


12.24 


67 


QSO+SB 


OM 


— 


— 


-1.42 


+0.48 


EL 


Li'pari ct al.2000b 


IRAS05024-1941 


-1150 


0.1920 


12.43 


— 


S2 


M 


— 


— 


-2.30 


-0.47 


EL 


Lipari et al.2000b 


IRAS05024-1941 


-1676 


0.1920 


12.43 


— 


S2 


M 


— 


— 


-2.30 


-0.47 


AL 


Rupke et al.2002 


IRAS05189-2524 


-849 


0.0420 


12.07 


— 


S2 


M 


— 


— 


-1.56 


+0.36 


AL 


Rupke et al.2002 


IRAS10378+1108 


-1517 


0.1360 


12.26 


— 


L 


M 


— 


— 


— 


— 


AL 




IRAS11119+3257 


-1300 


0.1885 


12.64 


— 


Sl+SB 


M 


1.12 


— 


-1.73 


+0.80 


EL 


Li'pari ct al.2000b 


IRAS11119+3257 


-950 


0.1885 


12.64 


— 


Sl+SB 


M 


1.12 


— 


-1.73 


+0.08 


OSM 


Zheng et al.2002 


IRAS13218+0552 


-1800 


0.2048 


12.63 


96 


QSO+SB 


OM 


— 


— 


-1.22 


+0.98 


EL 


Li'pari et al.2000b 


IRAS14394+5332 


-880 


0.1050 


12.04 


— 


S2 


MM 


— 


— 


-1.98 


-0.40 


EL 




IRAS15130-1958 


-780 


0.1093 


12.09 


— 


S2 


M 


— 


— 


-1.82 


-0.36 


EL 


" 


IRAS15462-0450 


-1000 


0.1001 


12.16 


— 


SI 


M 


— 


— 


-2.13 


-0.05 


EL 


" 


IRAS15462-0450 


-1110 


0.1010 


12.35 


— 


SI 


M 


1.32 


— 


-2.13 


-0.05 


OSM 


Zheng et al.2002 


IRAS19254-7245 


-800 


0.0597 


12.04 


30 


QSO+SB 


PM 


— 


— 


-1.70 


-0.11 


EL 


Colina et al.1991 


Mrk 231 


-1000 


0.0422 


12.53 


32 


QSO+SB 


M 


1.83 


Yes 


-1.49 


+0.11 


EL 


Li'pari et al.1994 


Mrk 231 


-1100 


0.0422 


12.53 


32 


QSO+SB 


M 


1.83 


Yes 


-1.49 


+0.11 


OSM 


Zheng et al.2002 


NGC 3079 


-1600 


0.0040 


10.49 


2 


L+SB 


s P 


— 


— 


-3.01 


-1.41 


EL 


Heckman et al.1990 


NGC 6240 


-930 


0.0243 


11.83 


15 


L+SB 


PM 


— 


— 


-2.16 


-0.40 


EL 




IR QSOoffset-OF 


























IRAS00275-2859 


-730 


0.2792 


12.64 


9 


QSO+SB 


PM 


1.47 


— 


-1.58 


-0.12 


OSM 


Zheng et al.2002 


IRAS00509+1225 


-1110 


0.0600 


11.87 


20 


QSO+SB 


PM 


1.40 


— 


-0.70 


-0.32 


OSM 


Lipari et al.2004d 


IRASZ01373+0604 


-1660 


0.3964 


12.30 


3 


QSO+SB 


- 


2.33 


— 


-0.50 


-1.59 


OSM 


" 


IRAS02054+0835 


-1355 


0.3450 


12.97 


— 


QSO 


- 


2.42 


— 


— 


— 


OSM 


Zheng et al.2002 


IRAS02065+4705 


-500 


0.1320 


12.27 


— 


QSO 


- 


1.65 


— 


-1.96 


— 


OSM 


" 


IRAS04415+1215 


-875 


0.0890 


12.41 


— 


QSO 


- 


1.74 


— 


— 


— 


OSM 


55 


IRAS04505-2958 


-1700 


0.2863 


12.55 


20 


QSO+SB 


PM 


1.33 


— 


-1.41 


-0.32 


OSM 


Li'pari ct al.2004d 


IRAS06269-0543 


-550 


0.1170 


12.49 


— 


QSO 


- 


0.57 


— 


-1.36 


+0.19 


OSM 


Zheng et al.2002 


IRAS07598+6508 


-2030 


0.1483 


12.45 


5 


QSO+SB 


M 


2.75 


Yes 


-1.32 


-0.04 


OSM 


" 


IRAS07598+6508 


-1920 


0.1483 


12.45 


5 


QSO+SB 


M 


2.75 


Yes 


-1.32 


-0.04 


OSM 


Li'pari et al.2004d 


IRAS09427+1929 


-1640 


0.2840 


12.61 


— 


QSO 


- 


2.98 


— 


— 


— 


OSM 


Zheng et al.2002 


IRAS10026+4347 


-650 


0.1780 


12.20 


— 


QSO 


- 


2.07 


— 


-1.26 


-0.88 


OSM 


" 


IRASZ11598-0112 


-970 


0.1510 


11.91 


— 


QSO 


- 


1.71 


— 


— 


— 


OSM 


" 


IRAS13305-1739 


-730 


0.1480 


12.21 


— 


S2 


OM 




— 


-1.24 


+0.21 


OSM 


Li'pari et al.2004d 


TRAGI QQA94-^Q^9 


Q9n 


n 1 7Qn 
u. J- ( yu 


12 49 
















ncu 

WOlVl 


7 Vibrio- of ol 9009 
Zjiieiig, et ai. zuuz 
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Notes. Col. 2: OF values obtained from the references 
included in Col. 9. Values between parentheses are possible 
detections. [O ill] MC indicates emission line [O in] A5007 
with multiple components. 

Col. 4: The Lir were obtained for [8-1000 /im], using 
the relation given by Sanders & Mirabel (1996). 

Col. 6: The properties of the nuclear activity were ob- 
tained from Veilleux et al. (1999), Canalizo & Stockton 
(2001) and Lfpari et al. (2003). SI: Seyfert 1, S2: Seyfert 
2, L: Liners and SB: starburst. 

Col. 7: For the morphological or interaction type we 
used the classification criteria of Veilleux, Kim & Sanders 
(2002b) and Sourace (1998); PM: pre-merger, M: merger, 
OM: old merger, MM: multiple merger, Sp: spirals. 

Col. 8: RF is the ratio of the emission line Fe 
HA4570/H/3 

Cols. 9 and 10: a are the IR colour indexes (for 3 IR 
bands). 

Col. 11: EL, AL and OSM indicate OF derived from 
emission, absorption lines, and offset emission line method 
(see the text), respectively. 
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Figure 2. Near-UV, optical and near-IR HST broad band images of the circumnuclear regions and the multiple shells, in Mrk 
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Figure 3. Near-UV, optical and near-IR HST broad band residuals images of the circumnuclear region and multiple shells, in Mrk 231 
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Figure 4. HST WFPC2 (I - B) colour images from the optical F814W and F439W filters. 
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Figure 5. WHT+INTEGRAL 2D optical spectrum of Mrk 231, for the Ha+[N n]A6848 blend. 
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Figure 6. WHT + INTEGRAL 2D spectra (for individual fibres) in the main structures of the nuclear and circumnuclcar regions of 
Mrk 231 (showing several OF components in the shells). The scales of flux are given in units of [erg X cm" 2 X s" 1 X A -1 ]. 
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Figure 6. Continued 
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Figure 7. WHT+INTEGRAL kinematics map of the ionized gas (Ha), for the central region of Mrk 231. 
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Figure 8. WHT+INTEGRAL maps of the emission line ratios for the central region of the main body of Mrk 231. The circle in panel 
(b) shows the nuclear region, of Mrk 231 (in the field of WHT+INTEGRAL), where the narrow component of emission lines [S n] are 
very weak. 
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Figure 8. Continued 
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Figure 9. KPNO ID optical spectra of Mrk 231. The scales of flux are given in units of [erg X cm 2 X s 1 X A x ] and wavelength in 

[A]. 
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Figure 10. Light curve variability of the Na ID BAL III system, between the years 1975 and 2000 (a), and between 1988-2000 (b). 
Showing also the best fit of the LC. For the fall of the LC the best fit is an exponential function. 
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Figure 11. UV ID HST FOS spectra of Mrk 231. 
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Figure 12. UV and optical spectra of the IR + GW/OF + Fe ii QSO IRAS 21219-1757 showing UV BAL systems. 
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Figure 13. UV and optical spectra of IR + GW/OF + Fe II QSOs (including specially the BAL IR QSOs). 
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Figure 14. Optical spectrum of Mrk 231 superposed to the spectrum of the SN 1979c (Branch et al. 1981). 
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Figure 15. IR colour— colour diagram for IR mcrgcrs/QSOs with galactic winds (from Table 4 of this paper) and for standard QSOs 
(from the PG sample; Boroson & Green 1992). 
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